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AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 



Length. 
Time__. 
Force__. 

Power. . 
Speed 



I 
t 

F 



P 
V 



Metric 



Unit 



meter 

second 

weight of 1 kilogram. 



horsepower (metric), 
fkilometers per hour. 
\meters per second. _ 



Abbrevia- 
tion 



m 

s 

kg 



k.p.h. 
m.p.s. 



English 



Unit 



foot (or mile) 

second (or hour)., 
weight of 1 pound 

horsepower 

miles per hour 

feet per second 



Abbrevia- 
tion 



ft. (or mi.) 
sec. (or hr.) 
lb. 



hp. 

m.p.h. 
f.p.s. 



2. GENERAL SYMBOLS 



W 9 Weight =mg 

g, Standard acceleration of gravity =9. 80665 
m/s 2 or 32.1740 ft./sec. 2 
w 

m, Mass= — 

I, Moment of inertia =mk 2 . (Indicate axis of 

radius of gyration k by proper subscript.) 
n f Coefficient of viscosity 



v, Kinematic viscosity 
p, Density (mass per unit volume) 
Standard density of dry air, 0.12497 



kg-nr 4 -s 2 at 



15° C. and 760 mm; or 0.002378 lb.-ft"* sec. 2 



Specific weight of 
0.07651 lb./cu. ft. 



'standard" air, 1.2255 kg/m 3 or 



3. AERODYNAMIC SYMBOLS 



s, 

Su>, 

o, 

b, 

c, 
b* 
S' 

v, 

2. 

L, 

D, 

D 0 , 

A, 

d t , 
c, 

R, 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 

Dynamic pressure =^pV 2 

Lift, absolute coefficient C L =^ 

Drag, absolute coefficient C D =^ 



D 0 



Profile drag, absolute coefficient Cd 0 =^ 



Induced drag, absolute coefficient C Di =^ 



Pj 

qS 



Parasite drag, absolute coefficient 0 D = 



Dp 
qS 



Cross-wind force, absolute coefficient C c = 
Resultant force 



C_ 
qS 



Q, 
a. 



VI 



C P , 



Angle of setting of wings (relative to thrust 
line) 

Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000; or for a model 
of 10 cm chord, 40 m.p.s., the corresponding 
number is 274,000) 

Center-of-pressure coefficient (ratio of distance 
of c.p. from leading edge to chord length) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute ^measured from zero- 
lift position) 
Flight-path angle 
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MECHANISM OF FLUTTER 
A THEORETICAL AND EXPERIMENTAL INVESTIGATION OF THE 

FLUTTER PROBLEM 

By Thkodokk Tiikodohskx and I. K. (Iahhick 



SUMMARY 

The results of the basic flutter theory originally devised 
in 1984 and published as N. A. C. A. Technical Report 
No. 496 are presented in a simpler and more com pU U 
form convenient for further studies. The paper attempts 
to facilitate the judgment of flutter problems by a syste- 
matic survey of the theoretical effects oj the various param- 
eters. A large number of experiments were conducted on 
cantilever wings, with and without ailerons, in the 
N. A. C. A. high-speed wind tunnel for the purpose of 
verifying the theory and to study its adaptability to three- 
dimensional problems. The experiments included studies 
on wing taper ratios, nacelles, attached floats , and external 
bracings. The essential effects in the transition to the 
three-dimensional problem have been established. Of 
particular interest is the existence of specific flutter modes 
as distinguished from ordinary vibration modes. On the 
basis of the concepts introduced, results that are apparently 
paradoxical could logically be brought into conformity 
with, the theory. In fact, it is shown that there exists a 
rather remarkable agreement between theoretical and 
experimental results. A simph method is prexeided for 
numerical calculations of the flutter speed by routine 
oj)erations, requiring no reference to the theory. Appli- 
cation is made to a complete numerical example. The 
matter of identifying possible types of flutter in an airplane 
and of determining the parameters is briefly discussed. 
A section treating the subject of forced vibrations of a 
wing in an air stream and the question of air damping 
in its relation to flutter is included. 

INTRODUCTION 

The theory of flutter. — The problem of flutter is 
passing through a period of rapid development. Full 
cognizance is taken of the value of the theory; a simple 
or an empirical understanding of this problem is not 
available and could, at best, be of value only to the 
investigator. An exact treat men! of the basic flutter 
problem in two-dimensional flow, involving the im- 
portant functions F and 0 relating to the air forces, 
was given by Theodorsen in 1934. (See reference 1.) 



These functions are simple combinations of Bessel 
functions; they have been rederived in related form 
by Cicala (reference 2) in 1935, by Kassner and 
Fingado (reference 3) in 193G, and also by Ktissner 
(reference 4) in 1936, who pointed out the identity of 
the functions. At about this time, Garrick (reference 5) 
also established a check on the general functions F and 
0 by comparing them with expressions by Wagner, 
Glauert, and von K&rm&n and Burgers for special cases. 

The system of equations as given hi the original 
paper is 

{A) aA+a^B+aC+pD+'^E+p^F+hG 



T l2 v, 



(B) all+a^l+ PJ+ p^K+ pL+hM+^K^C(k)Z=Q 

(C) aN+ ajO+ 0P+ j3^Q+ M+hS+ 2k^C(/c)Z= 0 
where A, B, C, etc., are given on page 10 of reference 1. 

z =r + b+(2- a h + ^ + ^ 



and (table 1) 
Putting 



where 



C(k)=F(k)+iG(k) 
a = a 0 e iiu>t ^ 



kv 



the determinant of the coefficients of a 0 , and h 0 
becomes 

Raa + ilaa Ra0 + Ha0 Rah + Hah 
Rba-{-lIba RbS-^-Ubft Rbh^Hbh 
\Rca + Uca RcP + ilcti R ch+ ll ch 

where the i?'s and I's are listed in the appendix. 
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This determinant put equal to zero contains two simul- 
taneous equations, the solution of which determines X 
and 1/k from which the (unknowns) flutter velocity and 
flutter frequency are obtained. Only the diagonal 
terms with bars contain the quantity X. All terms 
contain or are functions of 1/k. 

The work of numerous investigators, notably Becker 
and Foppl, has shown that the structural friction is 
mainly a function of amplitude, not of frequency. 
The structural friction can be described by a force in 
phase with the velocity but of a magnitude proportional 
to the restoring force. "With each res to ring-force 
term, say aC a , there will be a friction term ia g a C a , in 
which g a is the damping coefficient. The net result is 
very simply that the restoring-force terms aC a , $0$, 
and hC h have been replaced by terms of the form aC a 
(1 + ig a ) } pCp(l +igp), hCh(l+ig h ). These friction coef- 
ficients occur only in the diagonal I terms of the deter- 
minant. 

Technical flutter problems and the flutter param- 
eters. — Experimental evidence, some of w hich is pre- 
sented later, has been accumulated which indicates 
that, in the two-dimensional problem, the flutter speed 
can be closely predicted from the theory if the para- 
meters are given with accuracy. (In fact, it seems that 
in some cases the flutter speed can be used to determine 
some parameters more precisely than by a direct 
method.) In the two-dimensional problem of flexure- 
torsion-aileron flutter, about a dozen different quantities 
are required to calculate the flutter speed. The 
determination of these parameters requires technical 
skill and experience and is perhaps the most difficult 
step in the solution of the flutter problem. A knowledge 
of the functional dependency of the flutter speed on 
each parameter is essential in order to obtain sufficient 
accuracy in the determination of the important ones 
and to prevent waste of time on those of less influence. 
This need is partly the purpose of the material given in 
this paper. 

One of the problems in connection with an actual 
airplane is the identification of the combination of 
vibration modes that may cause flutter. In regard 
to wing flutter, in the case of flexure-torsion, the situa- 
tion is fairly clear. It will be shown that normally the 
most important parameter is the ecntcr-of-gravity 
location. This constant can be obtained with consider- 
able accuracy in the design stage. An accurate value 
of this parameter can also be experimentally obtained 
as the "dynamic" torsion axis, that is, the 4 axis around 
which the wing, owing to the low bending frequency, 
oscillates when put into torsional resonance. The 
location of the (static) torsional stiffness axis is much 
more difficult to calculate or to determine experimentally 
but fortunately, as will be observed, its effect on the 
flutter velocity is small. 

The internal damping coefficients are, moreover, of 
fairly small influence in flexure-torsion flutter; these 



parameters are also fairly difficult to obtain. On the 
w hole, however, it may be said that this case of 
Mutter can be fairly well handled. 

Another important case, that of the combination 
flexure-aileron, was shown by the original study 
(reference 1) to be an essentially different type of 
flutter from flexure- torsion. Its primary characteristic 
is that, normally, the Mutter is limited to a range of 
speeds. Below and above the extremes of this speed 
range there 4 is aerodynamic stability. A reduction of 
the static moment of the aileron with respect to the 
hinge (balancing) reduces the range; that is, the lower 
limit is raised and the upper limit is lowered. Damping 
in the structure is found to have the same general effect. 
Sufficient internal friction will, in fact, completely 
eliminate the danger of flutter, as will also complete 
mass balance or the proper combination of both. 
The structural friction of a wing system, although not 
readily predictable, can be obtained by a ground test. 

In regard to the tail assembly, the difficulty is some- 
what greater since it may not be easy to identify the 
most dangerous combination or to predict or even to 
measure the necessary parameters, including the 
structural damping. It therefore seems that ce.i haiti 
empirical or semiempirical aids will be required m<\ 
that it will be necessary, for a time at least, to resort 
to flight-test methods as a final assurance against tail 
flutter. 

The transition to the three-dimensional case of actual 
flutter is quite complex. It is necessary to consider an 
aerodynamic span effect (which fortunately is very 
small, see reference 6) ; the variation of the parameters 
along the span; the possibility of higher-order deflec- 
tion modes; and, in certain cases, fractional span effects, 
as for partial ailerons. The most promising manner 
of attack on such problems is by means of the two- 
dimensional treatment with the introduction of certain 
weight functions and average parameters in conjunc- 
tion with a study of representative models of reason- 
able simplicity, followed by a crystallization of the 
collected experience into generally applicable semi- 
empirical correction factors. The present paper makes 
initial studies with this purpose in mind. 

It is realized that, for high values of the flutter speed, 
a correction must be made for the effect of compressi- 
bility. In the first order, this effect is due to a change 
iu the slope of the lift curve. The air forces in the 
steady case are known to be increased approximately 
in the ratio M 2 , where Mis the Mach number. 

Consequently, a decrease in the flutter speed, roughly 
as (1— -M 2 ) 1/4 , is expected. This correction, although 
small through the usual flight range, becomes appreci- 
able for speeds near sound speed. Until experimental 
verification is available, such correction is preferable 
to none and should be applied for high-speed airplanes. 
(See footnote 2, p. 9, for details.) 



MKCJI.WJSM 

Content of paper. — A straightforward scheme is pre- 
sented in the first section for routine calculation of the 
flutter speeds in the two-dimensional types; case 1, 
flexure-torsion; case 2, flexure-aileron; and case 3, 
torsion-aileron ; and in the three-degrees-of-freedom 
type, flexure-torsion-aileron. A numerical example, 
referring to a modern largo airplane, is included. 

The second section deals with a survey of the effect 
of the flutter parameters on the critical velocity. 
The effect of changing the parameters within certain 
practical limits in cases 1, 2, and 3 is shown by a num- 
ber of charts. 

The discussion in the third section deals with the 
transition to a three-dimensional case, showing how a 
"ropi-esentative" two-dimensional wing may be used 
to give the essential results. Both uniform and tapered 
cantilever wings are included. The question of the 
probable occurrence of higher-order bending modes in 
(hitter is also discussed. The effects of "friction" and 
"coupling" are especially pronounced in higher-order 
flutter. 

It is pointed out that the deflection mode occurring 
in flutter is quite different from that of the static 
condition and that the lowest bending frequency 
involved in flutter is greater than that of the lowest 
ordinary vibration mode. A new concept of flutter, 
that the mode arising in flutter is such that the flutter 
speed is a minimum, is then introduced. In other 
words, if all primary variables including friction could 
be included in the analysis, the actual mode would he 
determined from all possible modes as the one giving 
the minimum critical speed. This concept is useful in 
explaining certain otherwise paradoxical, experimental 
facts. The extreme difficulty of a direct analytic 
attack on the general case, even if nil the physical 
parameters were specified, justifies the adapting of the 
two-dimensional treatment supplemented by empirical 
information obtained on actual w T ings. In fact, as will 
later be shown, the corrections are small. 

Almost 100 separate experiments were conducted 
in the 8-foot high-speed tunnel. The fourth section 
deals with the experimental tests and results. About 
one-half of these tests pertain to flutter of wings in 
flexure-torsion; the rest pertain to aileron flutter. 
Cantilever wings of aluminum and of built-up wood 
construction were used. The tests were performed on 

. a conveniently large scale, most of the wings having a 
chord of 1 foot and a span of about 7 feet. The air 

| speeds ranged from 50 to about 300 miles per hour. 
A number of safety devices had to be employed to 

I prevent the ruin of the tunnel equipment. 

A section is included showing the theoretical effects 
of the air damping on the forced vibrations of a two- 
dimensional wing system. This study leads to a more 
comprehensive understanding of the flutter condition, 
since it studies not only the critical speed but also the 
approach to this speed. A number of figures are pre- 

I sented that show the nature of the response curves in 
both one and two degrees of freedom. It is perhaps 



OF FLUTTER 3 

worth mentioning that von Schlippe (reference 7) has 
employed an experimental flight method for deter- 
mining the critical (hitter speed, which is based upon 
the use of an impressed alternating exciting force. 
The practical value of experiments of this nature is 
yet somewhat doubtful since the flutter usually comes 
on rather explosively. In any case, the theoretical 
results are of interest because they indicate the critical 
frequency as well as the growth of the maximum 
response as the critical speed is approached. 

METHOD FOR ROUTINE CALCULATION OF FLUTTER 

SPEED 

The calculation of the flutter speed can be reduced 
to a routine procedure by the following scheme. Noth- 
ing more involved arises than the calculation of the 
numerical values of double and triple determinants. 

Given are a maximum of seven original parameters 
k, r a 2 , a, x aj rp 2 , x$, c y from which are formed the 18 con- 
stants Aau A a2 , Ap u Aft, etc., defined as follows: 

J "=T+(s+" s ) 

^=£-^-(f-«>,] 

-4/33 = -'C^ + 7 7 lo) 
7T 

A hl = a 

K 

K 7T \ K 7T / 

B«=±(p-T,-iT,) 

7T 

0,2=1 

K 7T 

O»=0 

<?,„=-+ 1 
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These constants are obtained from the original varia- 
bles and from the T table (table 2) given at the end of 
the report. 

Another set of quantities R aa , I aa , etc., will be 
needed; their expressions are as follows: 

R aa = ~A al +(j- a 2 ) ^"(! + 

T ^ ( T> I -^ 12 -^ 10 2^ I Tl2 7n 0 



These six quantities are derived from the constants 
already given and from two additional quantities Fand 
G, which arc functions of l/k. The quantity l/k is, 
in reality, the independent variable in the problem. 
The quantities F and G occur in the forms 2F, 2F/k' 2 , 
and 2G/k; their values are given in table 1 for different 
values of l/k. In order to facilitate the calculation of 
these quantities, the parts depending on l/k are given 
in tables 3 and 4. 

Additional constants involved are the frequencies 
Q a , Qp, and Q h , defined under the scheme for each case, 
and three damping constants g a , g$, and g h . Generally 
all these constants are not simultaneously needed. 
The four cases will next be solved. 

Case 1. — The problem is given by two quadratic 
equations, for convenience referred to as the "real" 
and the "imaginary" equations. The coefficients of 
each are given in the calculation scheme presented in 
the following section. The coefficient of the first term 
in each equation involves the constants g a and g h , the 
coefficients of internal friction or structural damping, 
which are given as original constants. The coefficient 
of the second term of each equation involves the #'s 
and the .R's and I's, just defined. The constant term 
in each of the equations is obtained by the schematic 
arrangement shown in the calculation scheme; it is 
made up from certain constants A\ y B lf C x , and D x 
together with the quantities 2F,2F/k 2 , and 2G/k. The 
quantities A u B u C u and D x are simple determinants 
built up from the constants A aX , etc. 

The coefficients of the two equations must be calcu- 
lated for a fixed value of l/k; these coefficients are then 



substituted into the equations and the solution, that 
is, the value of X, is found. The real equation usually 
has two solutions, and the imaginary equation usually 
has one. The values of X, or preferably of -yjX, are 
then plotted against l/k, and the procedure is repeated 
until continuous curves representing the two equations 
are obtained. (Attained judgment or the knowledge 
of the solution of similar cases may considerably reduce 
the labor involved because it is then possible to choose 
reasonable values of l/k at the start. For wings and 
ailerons, l/k is usually less than 5, very often around 
1 or 2.) The point of intersection of the two curves 
represents the flutter point. Read off the values of 
X and I//-. The flutter speed is then given by the 
expression 

r a o)J) 1 1 
\k k V X 

Case 2. — The coefficient of the first term in each of 
the two quadratic equations again involves the con- 
stants of internal friction g fi and g h . The coefficient of 
the second terms is built up as in case L The constant 
term is built up likewise. Proceed as outlined for 
case 1. The critical speed is then 

<* h h 1 1 
Vk ^ -yj X 

where l/k and X are the values at the intersection point 
of the curves representing the real and the imaginary 
equations, respectively. There are usually two critical 
speeds. 

Case 3. — Case 3 requires a more laborious calculation 
of the constant terms; otherwise, the procedure is the 
same as for cases 1 and 2. The flutter speed is given by 

^ r a co a b 1 1 

-\/k k x 

Three degrees of freedom. The case of three degrees 
of freedom requires the solution of two third-degree 
equations in X. The constants of the first, the second, 
and the third terms are readily recognized as containing 
only quantities already used under cases 1, 2, and 3. 
The expressions for the constant terms of the two 
equations, D R and D 1 , involve three-row determinants 
but can be obtained by straightforward calculations 
for each value of l/k. The point or points of inter- 
section of the two curves representing the equations 
are again representative of the critical speed, which is 
given by 

_ r a a> a b 1 1 
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CALCULATION SCHEME 

Case 1 (h, a). — 
Real equation: 

Coefficient of A 2 : Q h tt a (l—g fl g a ) 

Coefficient of A": ^ h {Raa — (Jhha)+^ a {B ch —g a I ch ) 

Constant: M x *=A t +l^+c£p 

I magina ry equation: 1 

Coefficient of .V 2 : il h i2 a (g h + g a ) 

Coefficient of A r : %(R a <&h + Iaa)+toa(Rchga + Ich) 



Constant: M 7 =^( D x + d 



Ci=- 



A a i A hX 




C a \ Cm 




A a] - 


G+a) 


Gat 1 






1 


Chi 


Am A a 2 


Chi Ca2 1 


A a 2 A hi 




C a 2 Chi 





1 C hl 



A a2 -(i+a) 
C a2 1 



0, 



r a o) a b 1 1 



Case 2 — 
Real equation: 

Coefficient of X 1 : —<ltf}ii) 

Coefficient of X: QpiRch—gpIctd + tohiRbfi—gJbit) 

Constant: M 2 ^=^ 2 +Z 2 p+(5 2 +5 2 p^+C r 2 ^' 
Imaginary equation: 

Coefficient of A 2 : tyilhigp+g,,) 

Coefficient of X: ^fi(Reh9p+Ieh) + ^h(Rbfi0h + Ibfi) 

» Note that when the friction coefficients g arc zero, a factor \/k can be canceled out 
of all terras in all imaginary equations. 



Constant: ^i^j^+ft^-^+JStp)^] 



A>= 



B 0l B hl \ 

Cfil Chl\ 



B 2 = 



H-2 



a= 



A,= - 

T> T12 
t>fii — 



Chl 



2tt 
Cfii 1 

Ti2 
Cfiz 1 

1 



Z 7T 



^12 o 

2^ 
1 C M 



7 12 n 
~ £>hi 



IT 

I C h 



D 2 = 



Bgz B h i 

C/32 Cftl 



Ctfftft 1 1 

Case 3 (a, — 
Real equation: 

Coefficient of A 2 ; f2 a fy(l -# a ^) 

Coefficient of A r : Qa(Rbfi—gaIbfd+tofi(Raa—g?Iaa) 

Constant: M 3 ^=^3+^3p*+(^3+5 8 p)^ 

{< 

Iniii<riiiiti'y equation: 
Coefficient of A' 2 : il tt {lp{g tt -\-()tt) 
Coefficient of X: il. a (R, s g a +h fi ) + ik{R aa gg+I aa ) 
Constant: M/= ][z> 3 + Ap+ (^^p)^ 



^3 = 



-4al 


A* 




A a2 


Afi 2 


#al 


Bfiz 




B a2 
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B»=(i-a) 



(§+«) A* 
2* B » 



— (i+a) 
Bo1 2i 



£«=-(*-<*) 



C,= -(i-o) 



(| + a) ^ 3 
2t 



10 
7T 



A 2 -d+a) 
B " 2 2^ 



2* B * 



(i + a) A„ 2 
2tt 



7'n 



4. 



■(*+«) 
2tt 



-d+o) 
5 <" 2T 



•d + «) 4n 
2^ 



c 3 = 



•d+a) A& 

7>Z #03 



Aa2 Api 

B a2 Bpi 



A a i A ^2 
B a i Bp 2 

A a 2 Aft 
Ba2 Bp? 



Three degrees of freedom with friction. — 

Real equation: 
Coefficient of .Y 3 ; n a Q^ n (l—(j a gp—(j^j n —g h g a ) 

Coefficient of X 2 : fl tt fy,[(l — g a g^)R ch — (g a +9fi)Ich] + ^M(l— g^h)Raa—(gp+gh)Iaa] 
+Q h Q a [(l —g h g a )R b(i — (g h +g a )hfi] 

Coefficient of X: 0t(J£*-*JkW 

Constant: ^=«+^+(«+^)^+(r+^p^ 

I mamillary equation: 
Coefficient of X 3 : S^^^a+^+i/^+fa^/O 

Coefficient of X 2 : fi a fy[(i —g a g^)Ich+ (ga+gfi)Reh] + ^M(i —gfi9n)Iaa+ (ga+gh)Raa] 
+ V h Q a [(l—g h g a )I b(} + (g h +ga)Rbfi] 

Coefficient of X: Qa(M 2 I +g a M 2 *) + Q (t (M l I +g 0 M l *) + W 
Constant: Z^[tf+Z?^+(r+^^ 
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U= 



U= 



A a \ A$2 Am 

B a i Bq2 Bhl 

C a \ C$2 Cj,l 

A a2 Aft A hi 

B a 2 B)i\ 

C a 2 C/33 Chi 



A a2 A$i 

B a 2 Bpl 

( a2 (• (51 



Cm 



Note that five of the determinants occurring in the 
expressions for T and T occur also in the expressions 
for S and S. 

MSB 

r a co a b 1 1 

NUMERICAL EXAMPLE 



The following example refers to a modern large 
airplane. The parameters, which were furnished by 
the manufacturer, arc: 

r 2 =0.25 



k = 0.25 
a=-0A 
x a =0.2 



c = 0.6 
^'=0.0012 



A verbal description of the representative parameters 
used in the example is: wing density, <r, about 2.5 
pounds per square foot per chord length (in ft); 
stiffness-axis location, 30 percent of the chord from the 
leading edge; c en ter-of -gravity location, 40 percent of 
the chord from the leading edge; aileron length, one- 
fifth of the total chord; balanced aileron (center of 
gravity of aileron at hinge axis, xp = 0) . (The structural 
damping coefficients g a , g 0 . and g h will be kept zero, 
corresponding to a safety factor.) It is not necessary 
to specify the chord length 2b and the torsional fre- 
quency co a until the final step. The following fre- 
quency ratios, however, are specified: 



(2)'= 3/2 



That is, the torsional frequency is four times the bend- 
ing frequency and the aileron frequency is 1.22 times 
the bending frequency. The constants from which are 
composed all the determinants in the calculation 
scheme are tabulated as follows: 



A«\ = 1-285 
A*2= .9 
A 01= .02374 
A 02= .23679 
Afii= .40744 
Am-1.2 


B nl =0. 02374 
B«2= .04009 
Bfii= .007028 
B(i2= .021177 
Bfit= .01651 
B hi = .02322 


Cal = 1.2 
C„2=l 

Cffi= .02322 
C&2= .14238 
C<J3=0 
Cm = 5.0 




^- 2 =0. 00C357 

-7T 


1 



The equations are written explicitly for = 1, that 
is, 2G\k= -0.2006, 2F= 1.0788, 2^= 1.0788. 

Case 1 (flexure-torsion). — 

Real equation: 

The coefficient of X 2 is 

F6 X4X1 4 
The coefficient of X is 

i(- 1.285-0.12593) + 1 (-5.0 + 0.2006) = -5.15213 
The constants 

^=4.985 
^ = -0.125 

A = -3.3 

Hence the constant term is 

M 1 B =4.985+ (-0.125) (-0.2006) 
+0.7 (1.0788) =5.76524 

The real equation is then 

jx 2 -5.15213X+5.76524 = 0 
4 

Imaginary equation: 
The coefficient of Xis 



1 



(0.9 - 0.07703) + 1 (1 .0788) = 1 .28454 



The constant term M/ is 

-3.3 + 0.7(-0.2006)-(-0.125)(1.0788) = -3.30557 

The imaginary equation is then 

1.28454^-3.30557 = 0 

The roots of the real equation are X= 1.187 and I 
19.421, and the root of the imaginary equation is 
X=2.573, or 

-y]X= 1.089, 4.407, and 1.604 

These values of ^/Xnre plotted against 1/k in figure 1. 
The curves traced by plotting the roots are shown in the 
figure. The intersection is at ^X= 1.594, l/&=2.46. 
The flutter speed is then 

2.46 



^-2 X2X 1.594 



6co a = 1.5426co a 
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In the present example, the chord 2b is 12 feet and co a 
is 90 (corresponding to a torsional frc(picncy of Nf> ( .) 
cyles per minute); bu a is then 540 feet per second or 
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Figure h — Case 1. Numerical example. 
The roots -JXoi the real and the imagi- 
nary equat ions against 1/A:. 



O 12 3 

\/k 

Figure 2.— Case 2. Numerical 
example. The roots V X ot the 
real and the imaginary equa- 
tions against \/k. 



about 368 miles per hour. Hence the flutter speed is, 
for this case, 567 miles per hour. 2 

Case 2 (flexure aileron). — 

Real equation: 

The coefficient of X 2 is 



|X0.0012X1=0.0018 



The coefficient of X is 

0.0018(-5 + 0.2006) + l (-0.007028 + 0.020470) = 
0.004803 



2 The compressibility correction: Let the calculated flutter speed for the incom- 
pressible fluid be Vi and let the corresponding speed for the compressible fluid be v e 
Denote r,/c by Mi and v c /c by M e , where c is the velocity of sound. Then (see Intro- 
duction) 

1 



Mc 2 ~ 



Vl-Mc 2 



or, on solving for A/ r 2 . 



*M\ l ~—+—- ■ ) 



Forexample, with «?,=567m. p. h., Mi= 567/760 =0.746, Af f = 0.650, and r c =4 ( J4 m. p. h. 
Note that the example given refers to sea level; at altitude, the example should be 
based on another value of k and an appropriate value of the velocity of sound. 



The constants 

^2 = 0.034601 

^ 2 =-0.08255 

#2=0.008154 

£ 2 =-0.016510 

C 2 = -0.024980 

D 2 = -0.10258 
The constant term —0.07322. 
The real equation is then 

0.0018Z 2 +0.004803A^-0.07322 = 0 
Imaginary equation: 
The coefficient of X is 

0.0018X 1.0788 + 1(0.021177 + 0.0003184) = 0.023437 
The constant term M 2 r = —0.088554 
The imaginary equation is then 

0.023437-Y- 0.088554 = 0 

The roots of the real equation are X=5.182 and 
— 7,85 and the root of the imaginary equation is X= 
3.778 or (for the positive roots) V^=2.276 and 1.944. 
These values of ^JX are plotted against 1/k in figure 2. 
The curves traced by the roots are shown in the figure. 
Since no intersection exists, this case is stable. 

Case 3 (torsion-aileron). — 
Real equation: 

The coefficient of X 2 is 

1X/ 3 v/ 0.0012 nnnnyl . 
1 X 22 X -q-^- = 0 . 0 0 0 4 5 

The coefficient of X is 
l(-0.007028 + 0.020470) + 0.00045(- 1.285-0.12593) = 



The constants 
^3 = 0.008468 
^4 3 =-0.021110 
2*3=0.000799 
#3 = 0.002090 



0.012807 



C^= -0.003129 
C 3 = -0.004241 
A=— 0.026964 
Z> 3 = -0.001474 



The constant term M 3 /2 =— 0.021173 
The real equation is then 

0. 00045 X 2 +0.01 2807 vY- 0.021173 = 0 
Imaginary equation: 
The coefficient of X is 

(0.021 177 + 0.0003184) + 0.00045(0. 9-0.07703) 
= 0.021865 
The constant term M 3 7 = —0.030076 
The imaginary equation is then 

0~021865A r -0.030076 = 0 
The roots of the real equation are X= 1.567 and —30.03 
and the root of the imaginary equation is X= 1.375 
or (for the positive roots) <y]X=- 1.252 and 1.173. 
These values for the ^X are plotted against 1/& in 
figure 3. The curve traced by the roots is shown in the 
figure. Since no intersection exists, this case is also stable. 

Three degrees of freedom (flexure-torsion-aileron).— 
Real equation: 
The coefficient of X s is 

1X0.00045X^ = 0.0001125 
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The coefficient of X 2 is 

0.00045(-4.7994) + 0.0001125(- 1.41093) 
+#(0.Q13442)=0.001042 
The coefficient of X is 
l(-0.07322) + 0.00045(5.76524) + ^(-0.021173) 
= -0.07592 

The constants are 

#=-0.032848 jT=0.017042 
R= 0.077092 T=0.028790 
S= -0.004381 [7=0.103485 
~S= -0.000344 17=0.017720 

The constant term 7^=0.094635 

The real equation is then 

0.0001 125X 3 +0.001042A r2 -0.07592A^+0.094635 = 0 
Imaginary equation: 
The coefficient of X 2 is 

0.00045(1.0788) + 0.0001125(0.82297) + M(0.021495) 

= 0.005952 
The coefficient of X is 

l(-0.088554) + 0.00045(-3.30557) + #(—0.030076) 

= -0.097501 
The constant term D T = 0.1 1711 
The imaginary equation is then 

0.005952A r2 -0.09756lA r +0.1 171 = 0 

The positive roots of the real equation are X= 1.270 
and 21.0 and the roots of the imaginary equation are 
1.302 and 15.08, or -y/X= 1.126 and 4.58, and 1.141 
and 3.883. These values of ^X are plotted against 
1/k in figure 4. The curves traced by the roots against 
are shown in the figure. The intersection is at 
VA=1.06, 1/^=0.875. Hence 

1.06 

For 26 = 12 feet and co a =90, the flutter speed is 304 
miles per hour. 

These examples have been selected from several 
listed under the last part of the following section, to 
which the reader may refer for other examples, includ- 
ing the case of an unbalanced aileron. 

THEORETICAL SURVEY OF THE EFFECT OF THE 
FLUTTER PARAMETERS 

The purpose of tliis section is the study of the effect 
on the critical speed of the various independent varia- 
bles. Although the theory in itself permits the solution 
of any particular case without difficulty, it is somewhat 
difficult to obtain a perspective of the effects of the 
parameters. Because of the many variables, this survey 
has been limited to the magnitudes and the ranges of 
most practical interest. It is realized that the effect of 
increasing or decreasing a certain parameter is depend- 
ent on the values chosen for the others. As a mathe- 
matical experiment, it is possible to change one variable 
and to keep all the others constant. With reference to 
practical problems, however, the change of one param- 



eter is usually accompanied by unavoidable changes 
in several of the others. This fact must be kept in 
mind when actual or proposed changes intended to 
increase the flutter speed of airplanes are considered. 
This discussion is intended to give only the salient 
facts; the charts contain the complete data. 

CASE 1 (FLEXURE-TORSION) 

The flutter speed for case 1 is plotted in the coefficient 
form v/bco a . In the following graphs, the frequency 
ratio o) h /o) a is generally used as abscissa and the critical 
flutter coefficient v/bco a , as ordinate. 

The graphs under each of the following sections of 
case 1 are arranged in order of decreasing values of k, 





Figure 3— Case 3. Numerical ex- 
ample. The roots V X of the real 
and the imaginary equations 
against \Jk. 



Figure 4. — Three degrees of free- 
dom. Numerical example. The 
roots -y/x of the real and the 
imaginary equations against l/k. 



starting with k=1/2 (lightest wings) and ending with 
k=1/20 (heaviest). The range of k for present-day 
airplanes is approximately k=1/3 to /c= 1/15. The 
graphs are further arranged in order of increasing 
values of a, starting with the smallest values of a 
(stiffness axis in the most forward location). In most 
cases, the radius of gyration is kept at a fixed value 
r« 2 =l/4. 

Effect of center of gravity x a . — The effect of x a on 
flutter speed is given in graph I-A. It may be observed 
that there is usually a decrease in the critical speed as 
the frequency ratio o) h /o) a is increased from zero and 
that the curves tend to a minimum near the frequency 
ratio <a h /w a =l. There are cases, however, in which 
the minimum critical speed lies at w h /a) a =Q. The 
transition takes place for a certain small value of x a 
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(f) k=W, a = -0.2; v D /bo} a = l.l2. (g) k=M; a = -0.4; x a =0 stable; v D /b (l)a =2.2i. 

(h) k=M; a = -0.3; t;o/6wa=1.58. (1) k=M; a = -0.2; pd/6«««1.29. 

Graph I-A (a-i).— The effect of x a ; the flutter coefficient against the frequency ratio; r«2=i/4. Case 1 (ft, a). 
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This value is greater the larger the values of k (light 
wings). For instance, when /c= 1/4, a value of x a = 
about 0.1 (graph I- A (g)) brings the minimum near the 



speed near unity frequency ratio from zero to infinity. 
As may be observed later, structural damping will 
greatly alter the shape of the curve in this range. 




origin. For k i/lQ, x a must be close to zero (graph 
I— A (m)) to cause transition. The transition is critical ; 
graph, I A (m) shows that a 2.5-percent change in the 
position of the center of gravity changes the flutter 



The range of most practical importance is, however, 
the neighborhood of the zero frequency ratio. (For 
wings, the ratio is approximately 1/4.) In this range, 
the parameter of greatest significance is really the com- 
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bination %+a+x a . In other words, the flutter speed 
is very nearly a function of the location of the center of 
gravity with respect to the forward quarter-chord posi- 
tion and not of the distance relative to the stiffness 
Graph I-B (a) shows clearly that the value of a 



axis. 



light wings (jc= 1/3 to 1/5). Graph 1-C (d) gives very 
normal values of the parameters as used for most wings. 
These curves, for a given wing, may be taken to give the 
effect of altitude. Note that, for a given wing with 
k=1/5 at sea level, k becomes 1/10 at approximately 




6 .8 l.O 1.2 



Graph I-D (a-c).— The effect of structural friction; the flutter coefficient against the frequency ratio; *=Mo; a = — 0.2; j a =o.i. Case l (h, a). 



actually has no influence on the flutter speed. Outside 
of this range, that is, for larger values of co /4 /co a , the rela- 
tionship is less simple. Graph I-B (b) shows the 
dependency on the center-of-gra vit y locution for various 
positions of the stiffness axis a. For a constant x a , 
that is, for a constant distance between the stiffness 



15,000 feet, with a resulting increase in the Mutter 
speed under normal circumstances. For the case with 
X a =0.2 given in graph I-C (d), the increase in the 
flutter coefficient is from 1.(5 to 1.95, or about 20 percent. 
It is possible that, for very light wings, the Butter 
speed might decrease with altitude until a certain 




.2 



4 .6 .8 1.0 1.2 0 .2 A .6 

(a) (t->i. (b) *=Mo. 

Graph I-E (a-c).— The effect of radius of gyration; the flutter coefficient against the frequency ratio; a = - 



1.4 



(C) K=^0. 

•0.2, x a = 0.1. 



Case 1 (ft, a). 



axis and the center of gravity, the flutter speed is 
increased as the stiffness axis (and center of gravity) 
is moved forward. 

Graph I-C shows the flutter coefficient plotted against 
The normal range of wings is included in the 
diagram (the heaviest wings to the right). The dia- 
grams are arranged in order of increasing values of 
(c*>h/o> a ) 2 and of a. An interesting result is the existence 
of a minimum critical speed that falls in the range of 

207726—40 — 3 



altitude is reached. For high values of l/« (heavy 
wings), the flutter speed increases nearly as the square 
root of the wing density, 1/k. 

Effect of structural friction g a , g h . — Graph I-D is 
intended to show the effect of the structural friction on 
the critical speed. As the coefficients of friction are 
increased, there is a definite tendency for the often 
pronounced minimum flutter speed near av,/a> tt =1.0 to 
disappear and to produce response curves of the type 
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obtained for negative value of x a > In the range of most 
practical interest (w^/oj«^0), the torsional friction is 
the more important. 

Effect of radius of gyration r a . — Graph I-E is ar- 
ranged in conventional order. Note that the flutter 
coeflicienl in the low o h co Q range increases with in- 
crease in the radius of gyration. This increase in the 
ilutter coefficient does not necessarily correspond to an 
increase in the flutter speed; it does if the torsional 
frequency u a is kept constant. If the stiffness is 



(Values in the preceding table are given in relation to 
the value for k=1/5, r« 2 =l/5, which is the case of 
lowest wing density and smallest radius of gyration.) 
The speed corresponding to given stiffness drops if any 
mass is added so that r a , the density or both are 
increased. Hence, any mass added not for the purpose 
of increasing the stiffness or moving the center of 
gravity forward is detrimental. 

Flutter frequency. — The flutter frequency is shown in 
graph [-F. It is seen, for instance, that for small 




2.0 



(a) k=Yo\ a = -0A. 
(c) k=Mo; a = -0.4. 



(b) K=y$; a = -0.2. 
(d) k= Ho; a = -0.2. 



Graph I-F (a-d). — Flutter frequency ratio as dependent on x a against frequency ratio; r a - = Y\. Case 1 (h, «). 



kept constant, which means that co a is decreased as 
l/r«, the flutter speed is actually decreased, as is shown 
in the following table. 

FLUTTER SPEED FOR CONSTANT TORSIONAL 
STIFFNESS 

[a = -0.2, z«=0.1, (aW«a) 2 =0] 



K ^^^^ 


1/5 


1/4 


1/3 


1/5 


100 


97.6 


92.4 


1/10 


91.2 


88.6 


84.8 


1/20 


84 


83.2 


80.2 



values of co h /ca aj the flutter frequency is around 60 
percent of the torsional frequency co«; for higher values 
of the flexural frequency, the flutter frequency ap- 
proaches or exceeds the torsional. This graph is 
primarily of interest in connection with experimental 
flutter research. 

Coupling factor £. — Consider a two-dimensional case 
of flutter in which only a part of the total length of the 
(infinitely long) wing is given the second degree of 
freedom. This arrangement, because of the deficient 
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coupling, exhibits a higher critical speed. Call the 
fraction having both degrees of freedom, £. The results 
are shown for several values of £ in graph I-G. 

Divergence velocity and approximate flutter for- 
mula. — It can be shown that the divergence velocity 
may be expressed in nondimensional form as 



b0) a \ K l /2+a 

The divergence velocity r n !bo) a is given in graphs 
I-A. This velocity is usually higher than the Sutter 
velocity. 

An empirical expression, which is useful in quickly 
obtaining the order of magnitude of the flutter speed 
for small values of and which appears to hold very 
well for heavy wings (with k<1/10) is given by 



1L -It*. 



Graph I-B (a) shows the curve obtained from the 
empirical expression (dashed) and a curve based on the 
exact values (in full lines). 

CASE 2 (FLEXURE- AILERON) 

The flutter coefficient ,for case r ,2~is vjbw h . The fre- 
quency ratio a)p/co h is ordinarily used as abscissa. The 
graphs are again arranged in order of increasing wing 
density. Two values of the location of the aileron 
hinge axis c have been included. The first value, c= 
% } or the aileron chord equal to 25 percent of the total 
chord, is intended to represent a wing-aileron combina- 
tion; the second value, r = 0, or the aileron chord equal 



to 50 percent of the total chord represents a stabilizer- 
elevator or a fin-rudder combination. Several values 
of x$ and rp 2 and of the damping coefficients g$ and g h 
have been included. 

It should be mentioned that ordinarily, as shown in 
reference 1, case 2 differs basically from case 1 by the 
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.8 1.0 1.2 1.4 1.6 1.8 

Graph I-G.— The effect of the coupling factor {; the flutter coefficient against the 
frequency ratio; n=Mo;a = -0.2; x a =0.2. 

existence of a flutter range extending between a lower 
and an upper flutter speed. This range of flutter can 
be reduced or eliminated by various means. It is im- 
portant also to notice that, beyond a certain value of 
the frequency ratio co^/co,,, in fact, for a value slightly 
greater than unity, no critical speed exists, since the 
critical area does not extend much beyond this point. 
The reduction of the center-of-gravity distance from the 




.4 



.6 .8 



1.0 /£ 



1.4 1.6 



1.8 0 

CJ B /W h 

(a) Effect of ^ 

(c) Effect of 00 and gh', k=Ho; r<j 2 =Moo; xp = l Ao. 

Graph II-A (a-d).— Flutter coefficient against frequency ratio; c=y>. Case 2 (0, h). 



2.0 



(b) Effect of xp\ k= Ho; t&=\Uo. 
(d) Effect of rfc k= Ho; xp = ^o. 
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binge has the effect of reducing and finally eliminating 
the critical flutter area. Internal damping shows the 
same general effect. The fact that the aileron extends 
effectively over a shorter length is theoretically ex- 
pressed by a "coupling factor" £, which is the length 
of the aileron divided by the total length of the wing 
executing deflection. The effect of £ is shown in some 
of the graphs. 

Effect of frequency a)p/u h (c= l / 2 ). — Graphs 1I-A (a) 
and (b) show the effect of varying xp in reducing the 
critical area. The effect of damping is shown in graph 
II-A (c)and, finally, the effect of r} in graph II A (d). 

Effect of center of gravity xp (c= l / 2 ). — Graph II-B 
shows the flutter coefficient against the center-of- 
gravity distance xp ) giving, for two values of k, the 



be defined for each value of the frequency ratio. It is 
necessary then to choose the largest frequency ratio or 
the smallest unbalance, then to calculate the other 
value, and finally to choose the most practical combi- 
nation, using a margin of safety. 

Effect of radius of gyration rp (c=Y 2 ). — Graph II-C 
shows, for a typical w ing-aileron case, the effect of 
changing the radius of gyration for various values of 
the frequency ratio. 

Effect of frequency a^/co/, (c=0). — In the preceding 
graphs, the hinge axis was at c=%. Graphs 11 D. 
II-E, and II-F show the results for c=0. The curves 
are arranged in order and show the effect of Xp ) rp 2 , 
g h) and Qp for k=% and Xo. One curve is also included 
for k=% (graph II-D (d)). 




O .01 .02 .03 .04 

(a) c= x /z\K=\i; r/j 2 =M2o. 



.01 .02 .03 .04 



0 .01 .02 .03 .04 0 .01 .02 .03 .04 0 
(b) c=Vi\ K=\i\ r/j2=>6o. (c) c=^; k=Ho; r^=H2o. (d) c»?&; r^»Hoo. 



Graph II-B (a-d).— Flutter coefficient against x$ for various frequency ratios. Case 2 0?, h). 



effect of varying the frequency ratio o)p/<a h at three 
values of rp 2 . Note that for large Xp (beyond normal 
range) the type of flutter reverts to that of case 1 ; that 
is, the upper flutter speed becomes infinite for a certain 
value of xp. 

It is important to notice, by considering each curve 
in this figure, that Xp must be decreased below a certain 
value, which is rather critical, in order to avoid flutter. 
If xp is larger than this value, the lower flutter speed 
remains at a virtually constant, small value. The 
frequency ratio exhibits a similar effect; that is, flutter 
is eliminated beyond a certain frequency ratio often 
greater than unity, whereas for smaller ratios, the 
lower flutter speed remains at a low, nearly constant 
value. In other words, a critical frequency ratio can 
be defined for each value of the unbalance and, in- 
versely, a definite critical value of the unbalance can 



Effect of center of gravity Xp (c=0). — The figures 
are given in graph II-E, arranged as usual. 




.004 . 008 .012 . 0/6 . 020 . 024 . 028 . 032 . 036 



Graph II-C— Flutter coefficient against for various frequency ratios; c=l$; «=Mo; 
X0=Ho. Case 2 (p, h). 

Effect of coupling factor (c=0) — In graph II-F the 
effect of the coupling factor £ is shown for an extreme 
case of unbalance (xp large). The superimposed effect 
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0 .2 A .6 .8 1.0 IB 1.4 1.6 1.8 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 
(a) K=\i; r^=y 2 o; Xfi = Vio, V>o. (b) K=\i\ r^ 2 =^o; z/j=Mo. 



(c) k=4S; rj»«Ho; **=Mo, K>o. (d) r^=^o. 

Orapb II-D (a-d).— Flutter coefficient against frequency ratio, showing elTect. of x&, g$, and g h] c=0. Case 2 (0, h). 




Graph II-E (a-d).— Flutter coefficient against x& for various frequency ratios; c=0. Case 2 (0, h). 
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of damping is shown for the zero frequency ratio. 
Notice how the coupling factor (£-»0) gradually 
eliminates the flutter area, 

CASE 3 (TORSION-AILERON) 

Three graphs, III-A, III-B, and III C, are presented. 
There is a similarity to case 2. Graph III-A shows 
how the internal damping increases the lower flutter 
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0 .2 A .6 .8 1.0 1.2 1.4 1.6 1.8 

Qisph II-F.— Effect of coupling factor £; flutter coefficient against frequency ratio; 
c=0; k=H> rfi 2 =l&;ip=H (an extreme case of unbalance). Also efTect of friction for 

speed. Graph III-B represents data taken from an 
net unl ease of a light wing with a smaller aileron. 
Note the striking similarity to case 2. For the value 
^=0.0066 (completely unbalanced aileron), a)p/o) a 
must be greater than 0.6 to avoid flutter; for the more 
normal value ^=0.002, o^/co* need only be ^0.1. 
The flutter area is eliminated by reducing Xp to a 
slightly smaller value. 

Case 3 (torsion-aileron) is probably of less practical 
importance because the elimination ol* flutter for case 
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Graph III-A.— Effect of friction coefficients g n , gp; flutter coefficient against fre- 
quency ratio; c=0.5; n=yio', a = —0.4; fa 2 =H; xp=\io\ r/j 2 =Heo. Case 3. 

2 ordinarily excludes the possibility of flutter in case 
3; but it is noted that, in order to eliminate mass 
coupling in the torsion-aileron case, a complete balance 
of the aileron in the ordinary sense (^ = 0) is not quite 
sufficient. It is actually found in the ease of a heavy 
wing and no internal friction (with xp = 0) that the 
flutter speed is low, particularly near w ^=co a . Even 
a slight amount of friction, however, is sufficient to 
cancel the cause of this flutter. Graph III-C (fairly 
heavy wing) shows that, for no friction, a small over- 



balance (xp<C0) is necessary to eliminate flutter. 
For light wings, the effect is less pronounced and 
x 0 = O is usually sufficient. It may be observed from 
the original set of equations that true balance againsl 
rotation implies r^ 2J r(c— a) x$=0. 

THREE DEGREES OF FREEDOM 

In order to familiarize the reader with the complete 
case of three degrees of freedom and its relationship to 
the three subcases, a set of typical figures is shown. 
The constants used are the same as those in the numeri- 
cal example (p. 8) with some additions. Case 
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Graph lll-B. — Effect of ie; flutter coefficient against frequency ratio; c=0.6; *=H; 
a=-0.4; ra'-=>4; rs 2 =0.0012. Case 3(a, 0). 

1 is shown in figure 1 under the numerical example. 
The flutter coefficient v/b<a a =1.542. 

Case 2 is shown in figure 5; each part of the figure 
refers to different combinations of x$ and o^/cua. No 
flutter occurs for the combinations shown in figures 
5 (a) and 5 (b) because of the balanced aileron and 
none in figure 5 (d) because of the large aileron fre- 
quency. For the combination shown in figure 5 (c), 
there is a normal range of flutter with two flutter points 
shown. 

Case 3 is shown in figure 6; each part refers, respec- 
tively, to the same aileron parameters used in case 2. 
(Note that w A /w a is % in all cases.) The combinations 



5.0 




-.005 0 005 .010 .015 .020 .025 .030 



oc 0 

Graph III-C— Flutter coefficient against xp for various frequency ratios. c=0.5; 
k=Mo; ra 3 =M; r0 2 =Meo; a = -0.4. Case 3. 

shown in figures 6 (a) and 6 (b) are again stable because 
of the aileron mass balance. For the arrangement 
shown in figure 6 (d), the aileron frequency is not high 
enough to prevent flutter as it did hi case 2. Condi- 
tions are still worse for the combination shown in 
I figure G (c). 
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0 / Z30t230 /a30/23 

l/h 

(a) (^) , »P;%^0. m (j*y*m.Mp-*- (0) ^) 2 =0;. l<? =0.00(i6. (d) ( 5*) f -»<fc*» -0.0006, 



Ficure 5— Case 2. The roots -JX of the real and the imaginary equations against l/k. Same parameters as in numerical example except as indicated. 
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For the case of three degrees of freedom, figure 7 
shows the results arranged in the same order as under 
cases 2 and 3. For the conditions given in figures 
7 (a) and 7 (b), flutter existed only in case 1. The 
flutter point shown is therefore essentially case 1 
flutter. The value of the flutter coefficient, however, 



TRANSITION TO THREE-DIMENSIONAL FLUTTER 
PROBLEMS 

The previous theory relates to two-dimensional 
flutter and, strictly, to a wing of infinite length. The 
second restriction is not very troublesome, the aspect- 
ratio, or span, effect being relatively unimportant and 




has actually decreased from its case 1 value of 1.542 to 
0.70 and 0.825, respectively. 

For the arrangement shown in figure 7 (d), flutter 
exists in cases 1 and 3. Here the ranges completely 
merge, indicating stability at only very low speed. 
Flutter exists in all three cases for the combination 
8hown in figure 7 (c). The case 2 flutter can be rec- 
ognized, almost unchanged, while again the flutter 
ranges of case 1 and case 3 have merged, as in figure 
7 (d). 

Figure 8 has been included to show that there is a 
considerable lowering of the flutter speed for low values 
of the aileron frequency even though the aileron is 
balanced. This condition is probably not of primary 
concern because a small amount of friction, particularly 
g mt will restore the flutter speed to its full (case 1) value. 
It is to be noted, however, that a slight overbalance 
fe<0) may be desirable. 



by no means as great as the aspect-ratio effect asso- 
ciated with stationary flows. It may be disregarded 
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Figure 8.— Flutter coefficient against frequency ratio wpfaa. Three degrees of 
freedom. x=0.25; a = -0.4; z a =0.2; r« 2 =l/4; c=0.6; xp=0; r^ 2 =0.0012; g a =gK=0. 

and tacitly considered as a safety factor, since an air 
speed of the order of a few percent more than thai in 
two-dimensional flow is necessary to cause flutter. 
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Consider the case of a rectangular cantilever wing. 
Some authors have attempted a solution on the assump- 
tion that the response curves in torsion and deflection 
under normal conditions (zero air speed) may be used 
in the flutter theory. It is contended thai this assump- 
tion is false. Several rather interesting experimental 
results will be presented in the next section, which 
show directly and indirectly that the modes in flutter 
differ radically from the ordinary ones. The following 
questions arise: (1) Does the wing flutter in the first, 
second, or third, etc., bending "mode"? (2) Are these 
modes in any way related to the ordinary types of 
vibration modes? 

Consider first the ease of a very high bending mode. 
It is useful to consider the two-dimensional ease as 
representing t lu* "averages" of parameters and variables 
of the three-dimensional ease. The variable // now 
appearing in the (two-dimensional) equations refers 
to the "average" h which approaches zero even though 
the local It in the loops is very large. It is, furthermore, 
evident that the average curvatures are greater, the 
greater the mode. Both these conditions are equiva- 
lent to a Large coefficient of internal damping, since 
the work lost per cycle refers to a very small average h. 




AAA 

Figure 9.— Schematic figure for higher-order bending mode of cantilever wing 
restrained by wires and with deflections in phase. 

The coefficient, in fact, rapidly approaches an infinite 
value as the number n of the mode is increased. It is 
probable that the second "flutter" mode involves a 
coefficient of damping 20 to 50 times larger than that 
of the first. 

A study of the graphs with specific reference to the 
effect of damping shows that higher (hitter modes can 
he expected only under very special circumstances. 
This fact does not mean that flutter occurs in the low- 
est (zero air speed) bending mode. 

The bending frequency in flutter of a cantilever beam 
is determined by a certain minimum condition. The 
wing will, of course, flutter at the lowest speed possible. 
It will, therefore, not assume its lowest (stationary) 
bending mode but will tend to assume a mode of a 
higher frequency. Since this higher frequency tends 
to uncouple the h degree of freedom, the actual response 
ordinarily happens to be a cross between the first and 
the second modes. Large internal friction will tend to 
push the response closer to the first mode. The result 
is a flutter speed distinctly lower than that calculated 
on the basis of the frequency of the ordinary funda- 
mental bending mode. The flutter speed calculated 
by using the lowest bending frequency is too favorable. 
In the case of wings of small internal friction (solid 
metal wings), the act ual flutter speed is only about 0.9, 
the speed calculated on the basis of the lowest bending 



frequency. In the case of conventional wings, the 
error is apparently in the order of only 1 or 2 percent, a 
fortunate coincidence because it permits the use of a 
small experimental-empirical correction. This point 
of view is in harmony with the Rayleigh principle, 
which states thai any response function whatsoever 
corresponds to a frequency higher than that of the 
fundamental. 

To recapitulate: The bending frequency involved in 
the flutter of a cantilever wing is greater than that obserrt d 
at zero air speed; the more so, the lower the internal 
(lamping of the wing structure. 

This interesting phenomenon is demonstrated by the 
photograph of the flutter of a uniform cantilever beam 
(see fig. 18) presented in the experimental section, 
which shows that the maximum amplitude is not at 
the tip but is rather close to the middle of the span. 

Another very convincing experimental proof of this 
phenomenon, given in the experimental section of this 
paper, is that a counterweight at the tip section in 
front of the center of gravity actually lowers the 
flutter speed of a uniform cantilever wing. For a 
relatively small counterweight, the tip section is beyond 
a node in the h curve. In this same connection, 
another rather remarkable experiment was made: A 
cantilever wing flutters at about 200 miles per hour. 
The point wdiere the node of the second bending mode 
(at zero air speed) intersects the torsional stiffness 
axis was fixed by connecting this point by wires to the 
tunnel walls. The wing subsequently fluttered at 
150 miles per hour. The flutter stopped when the wire 
broke! The explanation is that the bracing wires 
"couple" a bending mode that was previously entirely 
"uncoupled." It should be noted that the frequency 
actually involved in this flutter is again in excess 
of that of the second bending mode (at zero air speed) ; 
large forces are therefore transmitted through the wire 
supports from the w r alls. 

In order to illustrate more convincingly how the sup- 
port wires lower the flutter speed, reference is made to 
figure 9, which shows a high-order bending mode of a 
wing. If this bending frequency is about equal to the 
torsional frequency, the low r est flutter speed is obtained. 
When the support wires are removed, the wing will tend 
to vibrate about a fixed mass center line, with the result 
that the average h deflection becomes zero and all h 
couplings disappear. The a moments and the h forces 
transmitted to the support are good measures of what 
may be called the effective values of a and h when the 
two-dimensional theory is applied to three-dimensional 
cases. For instance, the transmittal of a small h force 
to the support indicates that the positive and the nega- 
tive h values very nearly cancel. The h effect, although 
locally large, may very nearly cancel itself. This fact 
does not prevent the use of a certain (small) average or 
effective h in the calculations. With no internal fric- 
tion, the flutter speed is not changed. As was pointed 
out before, the use of the small effective h for higher 
modes is, in reality, equivalent to employing a greatly 
increased coefficient of internal friction. 
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Tli is discussion and associated experiments lead to the 
important conclusion: A bracing wire may lower the 
critical speed of a cantilever wing or fin. It usually does 
lower the critical speed when the internal damping in 
the structure is low. 

For a rectangular cantilever wing, there* is no difficulty 
in regard to the other parameters. With the bending 
modes considmul known, the variables a, j3, and li were 
given simply as average values and used in the two- 
dimensional solution. Inspection shows that the Butter 
speed of ;i uniform cantilever wing is essentially thai of 
the two-dimensional case involving the same parameters 
and the proper frequency ratio. 

A cantilever wing of normal tapered shape will nexl 
be considered. It is assumed that there is a similarity 
in construction along the span for each cross section. 
The mass is put equal t<> a constant times the square of 
tin 4 chord; static moments, to the third power: and the 
moments of inertia, to the fourth power of the chord. 
Further, the air force is proportional to the chord and 
the acting moments are proportional to the second 
power. 

Various weigh! factors of the form (6 b r )*f(x)dx are 
obtained, where / (./ ) is a weighted wing parameter and 
x is pleasured along the span. If the reference section 
is chosen in such a location that, for a particular n=m, 

which is always possible, then 

f(l) m+1 mdx^ff(x)dx 

In other words, the proper choice of a reference section 
renders ( he weight factors of approxima tely equal magni- 
tudes. If the reference section is taken too close to the 
tip of the wing, there will be a certain positive correc- 
tion; if chosen too close inboard, there will be a nega- 
tive correction. The correct value is thus virtually con- 
lined between definite limits. The most representative 
section will lie close to the three-quarter semispan 
location. 

In the two-dimensional case, the length along the span 
is considered to be equal to unity and this unity is 
treated as being large as far as span effects are con- 
cerned. If the length is different for the two variables 
considered, a slight modification of the theory is neces- 
sary. Kach length is considered to be long enough to 
permit disregard of aspect-ratio corrections for the air 
forces. 

This sort of consideration is of interest chiefly in the 
case of ailerons and tail surfaces. The equation giving 
the equilibrium of the ailerons refers only to the length 
of the aileron. The included area of the h curve is 
sometimes a small fraction of the total area under the 
h curve. This fraction will be called £ . 



The solution of the deflection-aileron case is given in 
reference 1 by 

_ \R bfi -f- il b p R bh + il bh I q 

\Bcfi-\-Hefi Rch + Hch\ 



and with the effect of £: 



Rbfi + UbP £(Rbh+iIbh) 
Rce+Ucp Rch+iLh 



= 0 



It is noticed that the factor £<0 describes a certain 
uncoupling of the system. The calculation of flutter 
speed can he performed for any coupling factor £. 
Again it should be remembered that the free-vibration 
modes are not identical with tin 4 flutter modes. A 
tendency exists for £ to approach unity since the 
aileron forces the motion of the wing. 

EXPERIMENTAL FLUTTER RESEARCH 

GENERAL 

The purpose of the experimental research was, first, 
to check the theory as regards accuracy and, second, to 
provide a basis for an understanding of problems met 
with in airplanes. 

These tests, about one hundred in all, were conducted 
in the N. A. C. A. 8-foot high-speed tunnel. (See 
fig. 10.) In order to protect the propeller, a heavy 
wire screen was inserted in the test section immediately 
behind the flutter model. For convenience, models 
having a flutter speed below 300 miles per hour were 
tested. 

The procedure followed was to increase the tunnel 
speed slowly until flutter appeared. If the flutter was 

of a violent type, the load was immediately dropped to 
save the model. In the tests on ailerons, the lower 
branch of the flutter curve was similarly obtained. The 
upper end of the range was obtained by the following 
method: The aileron was kept in place by restraining 
wires attached to its rear end and running across the 
tunnel. By manual operation of the wires from the 
outside, the arrangement could be conveyed through 
the dangerous ranges on slackening the wires, the 
operator would receive indication of incipient flutter 
until the speed had increased above the dangerous 
range. When the upper stable region had been 
reached, the wires were completely released and the 
conventional flutter-test procedure was reversed; that 
is, the tunnel speed was slowly decreased until the 
violent flutter appeared. The restraining win s were 
then immediately tightened, and the speed was noted. 
The effect of the very fine wire was shown to be negli- 
gible in the 1 released condition. 

DESCRIPTION OF WING FLUTTER MODELS 

All wings tested were cantilever wings and are based 
on the section given in the following table. 
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Figure 10 —Installation of wing 1 in 8-foot high-speed tun 
A MM OIL SHAPE USED IN FLUTTER TESTS 



Station % (per- 


Ordinate ?/ 
(maximum 


cent chord) 


thickness = 1) 


0 


0 


1.25 


.460 


2.5 


.600 


5 


.740 


7.5 


.820 


10 


.880 


15 


.940 


20 


.960 


f,0 


1.000 


40 


.960 


50 


.900 


60 


.780 


70 


.620 


80 


.440 


90 


.240 


95 


.132 


100 


.020 



The frequencies of the various wings are given in the 
main table of experimental data (table I). All section 
constants were obtained both by calculation and by di- 
rect testing. The basic section has its center of gra vity 
at 42.5 percent from the leading edge. The stiffness 
axis is at 32 percent but was artificially put at 30 per- 



nel. The stop shown was used in only a few initial tests. 

cent in the case 1 tests by chordwise cuts. (See figs. 
11 to 13.) 

In addition to obtaining the flutter speed of the plain 
wings, the effects of restraining wires, of mass balancing 
counterweights in various locations, and of large nacelles 
both at the wing and some distance away from it were 
studied. Experimental data are included in table I. 
In the aileron tests, the effects of mass balancing, hinge 
location, frequency, and friction were investigated. 

Wing 1 — Wing 1 (see fig. 10) was a rectangular can- 
tilever whig model of %-ineh duralumin plate of 12- 
inch chord by %-inch thickness by a free length of 6 
feet 9 inches perforated with closely drilled %-inch 
holes and covered by a %ooo-inch sheet of duralumin 
to give a smooth surface. The constants can be 
obtained from data in the experimental table. 

Wings 2, 3, and 4— Wings 2, 3, and 4 represent a 
series of cantilever wings of the same root section (1- 
foot chord by % inch thick), the same span (6 feet 9 
mches), bu1 having taper ratios, respectively, of 1:1, 
2:1 and 4:1. (See figs. 11, 12, and 13.) The wings are 
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made of duralumin and arc constructed to give simi- 
larity in strength and mass distribution. Note that 
the detail at or near the tips is a scaled-down replica 
of the detail at the root. The stiffness axis a is put at 
30 percent chord or a= — 0.4 by means of chordwise 
cuts. 

The three types of wing 2 (2A, 2B, and 2C) were so 
designated because the first one, 2A. finally showed a 




Figure 1 1.— Rectangular cantilever wing 2A. Note chordwise cuts used for purpose 
of lowering torsional frequency and for placing stiffness axis at 30 percent chord 
from leading edge. 

crack and had to be replaced with 2B, w hich is almost 
identical. Wing 2B finally broke al the root, was 
repaired by shortening it, and was used for some tests 
under the designation 2C. 

Wing 5. — Wing 5 was also a solid duralumin rec- 
tangular cantilever wing of 1-foot chord, 4-foot length, 
and 1-inch thickness at the maximum ordinate: it was 
used for aileron testing. (See fig. 14) Three ailerons 
were tested, 14, 24, and 34 inches long with 2, 3, and 
4 lunges, respectively. Most of the tests were per- 
formed on the 24-inch aileron (aileron A II). 

Tests were made for different spring-restraints on 
the hinge, with a balance counterweight on the out- 
board end (fig. 15) and with a special arrangement 
permitting the changing of the hinge axis from the 
forward edge of the aileron to about 30 percent of the 
aileron chord behind the center of gravity. 




Figure 12— Tapered cantilever wing Figure 13.— Tapered can- 

3; taper ratio both in chord and tilevcr wing 4; taper 

thickness is 2:1. Dimensional simi- ratio, 4:1. 
larity of cross section and cuts. 




Figure 14.— Cant ilever wing 5 used for aileron tests. 



MECHANISM 



OF FLUTTER 



Wings 6 and 7. — Wings 6 and 7 are model wings of 
normal density built by covering ;i balsa structure with 
Mo-inch mahogany. Wing 6 has the same external 
dimensions as wing 2 (fig. 16). Wing 7 has a root chord 
of 18 inches, a maximum thickness of 1.5 inches, and a 
taper ratio of 3:2 (fig. 16). All tapered wings were 
tapered equally in chord and thickness. 

DISCUSSION OF EXPERIMENTAL RESULTS 

The scheme already discussed of introducing flutter- 
bending modes completely fits the experimental results 
into the theoretical picture. Figure 17 shows the 
theoretical flutter speed for wings 2A, 2B, 3, and 4 
w ith the experimental points plotted. Wing 2A with 




Figure 15.— Wing 5 with aileron mass balanced by counterweight at outboard end. 

a flutter speed of 202 miles per hour obviously bends in 
a "first" flutter mode that approaches the second bend- 
ing mode in appearance and frequency (fig. 18). The 
flutter frequency calculated on the basis of this bending 
mode closely checks the measured (hitter frequency 
(fig. 19). Wing 3 checks equally well; its bending 
frequencies are noted in table I. Wing 4, the most 
tapered one, obviously collapsed (fig. 20) in the second 
flutter mode. (See fig. 17.) On this assumption, its 
experimental flutter speed also fits well in figure 17. 

Since the effect of the bending mode was brought so 
strongly into the picture, an independent study was 
made on the rectangular wing 2B and on the tapered 
wing 3 by attaching one point of the torsional axis 
rigidly to the tunnel walls by restraining wires. The 
results are shown in figures 21 and 22. Note that the 
wire attached to the tip had no effect on the flutter 



speed, which fact again tends to prove the contention 
that the flutter bending response is closely related to 




Fig UK e 16.— Rectangular wooden wing G (left) and tapered wooden wing 7 (right); 
taper ratio, 3:2. 

the second ordinary bending mode. Note also that the 
observed minimum speeds correspond very nearly to 



400 



300 



to 



200 



100 




.4 .6 .8 1.0 1.2 
Frequency ratio, Ctf^'/u^ 

Figure 17.— Theoretical flutter speed based on constants pertaining to wings 2A, 2B, 
3, and 4. k-Hq; a = -0A; x a =H; r<» 2 =0.3125; 6w a =38.8 miles per hour. Experi- 
mental test points are also shown, and flutter modes and frequencies are indicated. 



the minimum theoretical speed. (See fig. 17.) Of 
practical importance is the fact that a stay near the 
root of a wing gave a higher bending frequency and 



28 



REPORT NO. 685 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 




Figure 18 — Wing 2A in flutter, demonstrating first flutter mode. Note tendency for node at tip and maximum amplitude near middle. 



MECHANISM OF FLUTTER 



29 



definitely lowered the critical speed. There existed 
points, however, near the middle of both wings for 

1.6 



.6 





F/l 


Till 

ifter frequencit 
Theoretical 


i 1 




















O 


2A 
3 




















□ 

































































































































































































0 



.2 



A 



.8 



1.0 l.E /A 



1.6 



Figure 19— Theoretical flutter frequencies based on constants for wings 2, 3, and 4 
with experimentally observed values for wings 2A and 3. k = \W, a = — 0A;xa=\i\ 
r a 2=0.3125. Case l{h, a). 

which the stays caused the flutter speed to attain a 
large value. The explanation is that, with this point 



iixed, the average h value becomes very small and 
the h deflection becomes "ineffective." A relatively 
high flutter speed results. 

The matter of leading-edge counterweights has been 
investigated, in particular on wing 2C. Figure 23 
shows the effect of moving a counterweight along the 
span. The weight lias a rather surprising negative 
effect near the tip, indicating that, hi this case, there 
must be an h node inside the tip and again substan- 
tiating the theory of the flutter modes. Farther in 
along the wing there was an expected increase in the 
flutter speed. "When all three weights were applied 
at the same time, the flutter speed for wing 2C was 
increased to 2\)~> miles per hour, which IS in good 
agreement with the calculated value. 

A large nacelle at an inboard position (fig. 24) in- 
creased the flutter speed from 202 to 216 miles per hour 
w hen in the forward position and decreased it to 197 
w hen in the rearward position. 




Figure 20.— The effect of violent flutter (in second mode) on wing 4. 
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Large bodies disposed at some distance from the 
wing, such as floats, were very detrimental as regards 
t be critical speed. (See fig. 25.) 

Wing 1 showed a flutter speed that is in agreement 
with the predicted value within about 1 percent. This 
agreement is due to the considerable internal damping 
of tins wing. Wing 6, a rectangular whig of the same 
plan dimensions as wings 1, 2, 3, and 4, but of low 
density, showed a flutter speed about 3 percent below 
the theoretical value based on the measured parameters 
and the lowest ordinary bending mode. This result 
indicates that, for damped, low-density structures, the 
flutter mode approaches the first bending mode some- 
what more than hitherto indicated. 

Wing-aileron flutter lias been studied on wing 5. 
(See fig. 14 and table IA.) The theoretical response is 
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Figure 21 — Experimentally observed flutter speed as depending on location of 
restraining wire along axis of wing 2B. 

shown with proper constants hi figure 26 for the most 
representative aileron All, upon which most of the 
tests were made. A number of test points have been 
directly plotted in tin's figure. In order to obtain 
internal friction, a lead hinge was used in some tests. 
It is rather remarkable how well the theory is re- 
aflirmed by the test data. Apparently, if all param- 
eters could be satisfactorily determined, no flutter 
testing would be necessary. Tests in winch the hinge 
axis (fig. 27) was changed show the beneficial effect of 
decreasing the aerodynamic moment around the hinge. 
The lower flutter speed, which is the one of practical 
importance, is considerably increased as the hinge 
axis is moved backward. This increase is not only a 
center-of-gravity effect but is also caused by the de- 
crease in the aerodynamic moment around the hinge. 
Note that, as the center of gravity is approached, the 
flutter speed rather suddenly becomes infinite. 

AIR DAMPING OF FORCED VIBRATIONS 

This report has heretofore been concerned with 
a study of a border velocity separating stable and 



unstable velocity regions. Further light on the whole 
matter of flutter is given by a study of the vibration 
response of the wing system to impressed forces and 
moments, that is, generalizing the point of view 
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Figure 22.— Experimentally observed flutter speed as depending on location of 
restraining wire along axis of wing 3. 

from free vibrations to forced vibrations. Instead 
of the homogeneous system of equations (A), (B), 
and (C) (see Introduction), impressed exciting forces 
and moments introduced on the right-hand side of 
these equations are considered. In equation (A) a 
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Figure 23 —Effect of counterweights. Wing 2C. 



80 



term M a e i(wWo) /M6 2 , in equation (B) a term 
Mpetw+w/Mb 2 , and, in equation (C) a term 
PtfiW+**)/Mb are introduced. Here M a and Mp 
are the magnitudes of the sinusoidal impressed torques 
in the a. and 0 degrees of freedom, P {) is the magnitude 
of the impressed force in the h degree of freedom, co is 
the circular frequency of the forced vibrations, and the 
^'s are certain phase angles. 
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Figure 26.— Theoretical flutter speed bnsed on constants pertaining to wing 5 with aileron AIL (C=H; ^=0.0076; r^ 2 =0.0019; k = Mos; 6o>a = 30.5 in. p. h.) Experimental 
values are shown; flutter region is shaded. Dashed (theoretical) curve corresponds to friction coefficient ^ = 0.0125. 



It /t {) 



CASE 1 (h, a) 

By substitution of a=a 0 e^ t+ ^\ e««'+*>, in 

equations (A) and (C) and solving for (j^ ?V2 and 
a 0 e iipo (put M a =Pohm), there results 



25 percent of the chord 
m=-0.5. 



in front of the elastic axis, 



Raa + U, 

P 0 e^ Rca+iL 



aa me 
1 



f(^ 0 — ^2)1 



a 0 e^= 



in 



A^+iA/ 
meWo-^Ran+ilat, 1 

1 Tt eH + iIeh\ 



"6Mcco 2 V + iA/ 



from which both the amplitudes and the phases may be 
obtained. The R's and the /'s are listed in the appen- 
dix and Ai R and A/ represent the real and the imaginary 
equations listed under the calculation scheme for case 1. 

Consider the equation for h^e*** and denote the 
determinant in the numerator by N b) i. e., 



Nn= 



Raa + Ha. 
Rca + Ucc 



1 



If the excitation is only in the h degree of freedom 
m=0, i. e., there is no impressed torque about the 
elastic axis. If a single exciter were placed, for example, 
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Figure 27.— Flutter speed as depending on position of aileron hinge. Observe that 
no flutter exists when the hinge approaches the center-of-gravity location. Aileron 
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Figure 28— Photograph of model showing forced vihration response in two degrees 
of freedom as depending on air speed. Note flutter point. 

It is convenient to define a certain static deflection 
h st , which is the deflection due to force P 0 , 

Po Po 



Define a st as the static torsional response to an impressed 
moment P 0 b 



h st = 



C h Ma> h 2 



_Pob 



Then 
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Figure 28 is a photograph of a three-dimensional 
model of the response ratio h 0 /h st as a function of the 
exciting frequency ratio w/a> a and of the speed ratio 
v/bo)a for a case of deflection-torsion for which a critical 
flutter speed exists. In this example, the exciting force 
has been assumed to act in the deflection degree of 
freedom alone (m=0). At zero speed there exist two 
resonant frequencies corresponding to the natural 
frequencies in the separate degree of freedom. The 
air damping due to speed is zero and the response is 
infinite at these frequencies. (With friction g a and 
g h) the responses are approximately l/g Q and l/g h , 
respectively.) As the speed increases, the air damping 
increases and the response diminishes until, along one 
frequency branch (the a> a branch) the response becomes 
negligibly small. Along the other frequency branch, 
however, a minimum response is reached, after which 
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Figure 29.— Peakj-esponse ridges of figure 28. 
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A = V(V) 2 +(V) 2 



This result gives the steady-state deflection response 
h 0 in terms of the static deflection h st due to an impressed 
force or moment. The results of some numerical 
calculations will shortly be presented. 

The torsional response can be similarly calculated. 
Let N a represent the determinant in the numerator of 
the expression for a o e ,>0 : 
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Figure 30.— Tlic roots of the real and the imaginary equations (pertaining to case 
shown in fig. 28) against Ilk. 



34 



REPORT NO. 685 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



the air damping decreases and the response increases 
(usually rapidly) until the flutter point is readied. The 
critical flutter frequency lies between the two critical 



example cited, is shown in figure 30. The value of 
A/ is then obtained for various selected points 
(I /k, o>/co«) along Ai*=0. These values may be em- 
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Figure 31— Peak response ho/h tt against velocity v/bua for the case k= Mo; a=— 0.2; x a =0.1; m=0; («*/««)' ~H 8 • 

dashed curve, to gh^g a =0.l. Case 1 (h, a). 



1.6 



J. 8 2.0 



The full curve refers to friction coefficients gh = g a =0; the 



frequencies that exist at zero speed. The two peak 
response ridges for this case are shown in figure 29. 
These curves, it may be observed, illustrate the essential 



ployed to determine A. The numerator is easily 
evaluated for these same values of 1/k, co/co«. This 
process determines the peak response with sufficient 




.6 .8 1.0 1.2 1.4 
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Figure 32.— Peak-response ridges for three values of friction coefficients g a = gh=Q, 0.05, and 0.10. {k = \[q; a = —0.2;x a =0A\ («A/coa) 2 =K>)- Case I. 



characteristics of the three-dimensional figure. The 
procedure of calculation is as follows: 

The equations A^— 0, A/=0 are solved for co/oo a for 
various fixed values of 1/k, as already discussed, in 
order to locate the flutter point. This plot, for the 



accuracy for a given value of 1/k and conversion 
to speed is obtained by the relation v/bu a = (co/a> a ) 
(1/k). Such response curves, calculated with and 
without the effect of friction, are given in figures 31 
and 32. 
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Figure 33— Forced vibration response for deflection 
degree of freedom (gh=0). 
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Figure 34.— Peak 
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One degree of freedom (deflection). — Further light 
on the two-degree-of-freedom case, deflection-torsion, 
may be obtained by a discussion of the one-degree-of- 
freedom cases, deflection alone and torsion alone. 

When h—h 0 ei( O)t +* 2) is substituted in the deflection 
equation of motion, 

On the introduction of h st =P 0 /C h (static response in 
deflection to impressed force P 0 ), it follows that 



h st K\U) / 



1 



(R c n 2 +IJ)> {(T h -A h y+ (g h i\+B h yy 



response hojh.t against velocity ratio v/bu h for two values of the friction coefficient 
gh=0 and 0a=O.1. n=\io. 

The maximum response is then 



\hst/u 



B h +g h A h 



where 



A-l+«(l+f) 



2F 

g h is the friction coefficient. 

It is observed that the speed is determined by v/ba h = 
(«/«*)(!/*)■ 

The resonance response is obtained by putting 
b/bT h \ho/hst\=0 and solving for r». There results 



or 



A h *+ B h * 
A h —g h B h 

/oA f A h -g h B h y 
\A h *+B h 2 J 



Figures 33 and 34 pertain to this case. The results 
may be summarized as follows: The resonant frequency 
is practically constant and is approximately o) = cj h . 
The air damping at the maximum response is propor- 
2F 

tional to B h =K -y or, since the frequency is nearly 

constant, the air damping at the maximum response 
is proportional to kv. Away from the resonant fre- 
quency, however, the response quickly becomes inde- 
pendent of k. No flutter or self-excitation exists in 
this case. At zero air speed, the maximum response 
ratio is simply l/g h . 

One degree of freedom (torsion). — When a= a 0 e i(o}t+tpo) 
is substituted into the a equation of motion 

a 0 eHRaa+iIaaW= jjjfe e*** 

On the introduction of a st =M a /C a (static response in 
torsion to impressed torque M a ), it follows that 

1 



where 



r 

B 



[(r a -A a y+(g a r a +B a y]i 



=i+ 



g tt is the friction coefficient. 
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Figure 35. — Forced vibration response for torsion degree of freedom (axis at a=0; 
Klr a *=lio;g a =0). 

It is observed that the speed is determined by 
v/bG> a =(a>/a) a ) (1/fc). Put &/dr a |ao/a, t | = 0. Then 



r res =(°^\ 2 - A °+ B ° 2 

\ CO / res A a — g a B a 

AcA _/ A a -g a B a y 

W rM \A a * + B a *J 
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Figure 36.— Peak response ao/a,< against velocity ratio v>lbu a for two values of the 
friction coefficient g a = 0 and 0 a =O.l. a=0; »c/r a 2 =^o. 

of the torsional axis. The air damping at resonance is 
essentially proportional to k//v and to the wave- 
length parameter l/k. For the quarter-chord position 
of the axis, a = —0.5, the response a Q /a st is very similar 
to the deflection response h 0 lh st in the preceding case 
(fig. 33). For any position of the torsional axis back 
of the quarter-chord point, however, a peculiar result 
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Figure 37.— Peak response ao/a»t against velocity ratio v/ba) a for two values of a (0 and • 



7 8 9 /0 

}4) and two values of k (Mo and J^oo). Case 1 (h, a). 



The maximum response is then 

(A a 2 +B a r> 

+g a A a 



Figures 35, 36, and 37 pertain to this case. The 
resonant frequency is strongly affected by the position 



is obtained. The air damping increases with increase 
in speed and the resonant frequency decreases; but a 
speed is ultimately reached where the response increases 
again until, at a vanishingly small resonant frequency, 
the response is very large. In figure 35, which illus- 
trates a case for the midchord position of the axis, 
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Figure 38.— Peak response holh,, against velocity ratio v/buh for a case of flexure-aileron flutter, showing 
effect of changing friction coefficients gh=ga=Q to gh=gp=Q.2. (Other parameters are c=0.5; k=Ho: 
r/3 2 =Meo;^=Ho; («0/wjO 2as ^). 
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Figure 39.— Forced vibration response of aileron alone 
against frequency and velocity; friction coefficient (7/3=0. 



a=0 y the peak response occurs at v/bco a =7 .10. This 
sort of instability has been called divergence. The 
divergence velocity can be easily calculated as follows: 
Let w/w a — >0, >oo in the expression for a 0 la st . 



hen 



Limit / a \ 1 



or 



(JL) 2 - r ° 2 l l 2 



+a 



CASE 2 03, ft) 

A study of the response characteristics to forced vi- 
brations is also of some interest in the ease of flexure- 
aileron. The details are omitted here. Two figures 
are presented. An illustration showing the peak re- 
sponse ratio fi 0 /h st in this case with and without friction 
is presented in figure 38. A response for one degree of 
freedom of the aileron alone is shown in figure 39. 

REMARKS ON FLUTTER IN AIRPLANES 

WING FLUTTER 

The wing may flutter as a whole in torsion-flexure. 
This case is the most easily treated Experience with 
models indicates that this flutter speed may be calcu- 
lated on the basis of the measured constants with an 
accuracy of a few percent. The actual bending fre- 
quency involved in flutter is apparently not exactly the 
lowest ordinary bending frequency but a slightly 
higher value. 

Probably the most common type of wing flutter is 
case 2 (flexure-aileron). This type, as well as that 



involving torsion-flexure, is evidently symmetrical with 
respect to the fuselage. The ailerons would therefore 
be in phase and have a frequency considerably in excess 
of the wing-bending frequency. This condition is 
favorable. Any slack in the aileron cables, however, 
permits a motion that may cause a mild type of flutter, 
which should not be permitted for too long a time. 

A nonsymmetrical aileron motion would involve a 
second bending mode (nonsymmetrical). It is prob- 
able that, in most cases, the node would be close to the 
middle of the aileron and therefore poorly coupled. 

There remains to consider a complete case of flutter 
(torsion-flexure-aileron). Apparently cases do exist in 
which this type would appear at the lowest speed. The 
effect of the additional degree of freedom can probably 
be taken care of by a safety factor applied to the flutter 
speed obtained for two degrees of freedom. The cal- 
culation of the case of three degrees of freedom is per- 
fectly straightforward although more lengthy than the 
simple cases. 

TAIL FLUTTER 

In regard to tail flutter, the situation is more complex. 
The possible combinations are subdivided as follows 
into three main groups, which will be separately 
analyzed: 

(1) Vertical flexure of tail assembly. 

(2) Horizontal flexure of tail assembly. 

(3) Torsion of tail assembly. 

Vertical flexure. — It is possible, in general, to iden- 
tify two responses in vertical flexure ; one corresponding 
to the fundamental bending mode of the fuselage and 
the other, to the bending mode of the horizontal-fin 
arrangement. The frequency of the fin arrangement is 
slightly greater than the one obtained with the rear end 
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of the fuselage fixed in space. It is probable that the 
fuselage bending mode need not be considered. In any 
case, the flutter spree I calculated for each of the two cases 
will not differ very much because the density involved 
is about in the same ratio as the squares of the fre- 
quencies involved. Only the vertical fundamental 
bending frequency of the horizontal fin will therefore be 
considered. This bending frequency may couple with 
fin torsion and elevator motion. This motion is neces- 
sarily symmetrica] and simulates the motion of the 
main wing system. Since the elevator has no particular 
restraint to this motion, it is evident that an unbalanced 
elevator is highly undesirable. As in the case of the 
wing system, the most frequent cause of flutter is also 
the flexure-aileron combination. 

Horizontal flexure. — Horizontal flexure affects the 
vertical fin or fins and may be separately considered as 
a cantilever wing with an aileron. The rudder, which 
takes the place of the aileron, has no particular restraint 
to this motion and must therefore be carefully mass 
balanced . 

Torsion. — Torsion is composed of a relatively low- 
frequency type of flutter involving the fuselage and a 
higher frequency type involving the fin tips, which may 
be considered as fin flexure in opposite phase. This 
type of flutter is not common because the great stiffness 
of the torque tube prevents its occurrence (ft large). 
The rudder is a fleeted in the same manner as for 
horizontal flexure. 

In summary, it may be said that the expected causes 
of flutter are the in-phase motion of the horizontal fins 
as flexure-elevator and the motion of the vertical fin in 
flexure with rudder motions. 

With two rudders disposed at the ends of the hori- 
zontal fins, care must be taken that the flexure fre- 
quency is sufficiently high. The mass of the rudders 
at the ends of the horizontal fins also affects the param- 
eters used in the fin-stability calculation; that is, the 
density of the fins and the radius of gyration are in- 
creased. The center-of-gravity location may also be 
changed. 

GROUND TESTS 

DETERMINATION OF CENTER-OF-GRA VITY LOCATION 

From the theory, it may be observed that the location 
of the torsional stiffness axis is of fairly secondary im- 
portance. The location of the center-of-gravity axis, on 
the other hand, is of great importance. The application 
of a very low-frequency (zero) torque will rotate the 
wing around the torsional stiffness axis a; the application 
of a very high frequency (infinity) torque will cause 
the mass center line to remain stationary. As the tor- 
sional frequency for wings is several times larger than 
the lowest bending frequency, it can be shown with all 
desired accuracy that the axis observed for the tor- 
sional frequency is the center-of-gravity line. 



Assume the wing to be vibrating around an axis at a 
distance d behind the stiffness axis. The moment of 
inertia reduced to the center of gravity is 

I=M(r 2 -x 2 )b 2 
\ loving the axis to d increases the moment of inertia to 

M(r/-j a 2 )6 2 +M(x a -rf) 2 6 2 =M(r a 2 +(/ 2 -2r/^6 2 
The corresponding torsional stiffness is C a -{-d 2 C h . 
The frequency is consequently 

2 _ (C a ±d 2 C„W 
03 ~M(r a 2 +d 2 -2dxJb 2 

Then 

2 log co=log (C a +d 2 C h )b 2 -log Mb 2 (r a 2 +d 2 -2dx a ) 

The whig will assume the state of vibration giving the 
greatest frequency. By derivation, with respect to d, 

dCf t d x a p. 

C^+Mf h ~ rj +d 2 - 2dr„ ~ 

or 

d—Xoc 

Then 

d= *"i -(«»/«)* 

or, with d known, 

x a =d[l — (o)„/o)) 2 ]^d for co/j/co — >0 

That is, the center-of-gravity axis is slightly ahead of 
the dynamic axis (assuming both axes to be normally 
behind the stiffness axis). If the torsional frequency is 
very large, they coincide (d=x a ). (If the torsional fre- 
quency is very low, d=0, giving the stiffness axis.) 

In other words, the center-of-gravity location along a 
finished wing can be determined by establishing the 
dynamic torsional axis. 

DIMENSIONAL CONSIDERATIONS 

Proportionally increasing all dimensions of a wing 
while retaining all details lowers the frequencies in in- 
verse proportion to the size. The reference speed a>6 
therefore always remains the same, as do all other pa- 
ra meters including the wing density. The actual flutter 
speed therefore depends on the shape but not on the size. 
It is important to keep in mind, however, that the 
reference is to wings or tails similar in all respects. In 
reality, a lighter contruction is necessarily employed in 
larger wing sizes, resulting in a weaker structure and a 
general lowering of the critical flutter speed. 

The foregoing considerations are significant in the 
testing of models. Thus a true model constructed of 
the same material as the full-scale airplane will have 
the same flutter speed. For testing purposes, it is 
very desirable to have a fairly low flutter speed. This 
end may be achieved by employing models of special 
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null trials related to celluloid, which have a value of 
-y/E/p nearly five times smaller than that of materials 
normally used in airplane construction. 

The most desirable condition would he to use a 
material with the same density p as the airplane and 
with the moduli E and 0 } say, l/n times the original 
valued. Fortunately, the density of the model wing 
can he very simply corrected by using a suitable thick- 
ness of the materials. Thus, if the density of the 
material used is three times lower than that of the 
original, the thickness of the skin and all the internal 
members is increased by a factor of 3. 

It. should he further noted that the model can he 
critically checked as to accuracy of reproduction by 
direct measurements of its mechanical properties. In 
other words, all the parameters, including the reference 
quantity <ab, are directly measured on the model itself. 
The value of ub is usually close to the predicted value. 
The important point, however, is that it is not neces- 
sary to depend on a predicted theoretical value. 

Thus the feasibility of conducting direct flutter tests 
on models of actual airplanes or of its component parts 
is indicated. Some work of this nature is now being 
undertaken. The procedure may be of value in cases 
that are difficult to treat theoretically and should be 
of value in accumulating useful experience on special 
designs. 

GENERAL CONCLUSIONS 

L. The two-dimensional theory has been verified 
within the limits of error in the determination of the 
primary parameters. 

2. The most essential three-dimensional effect is the 
occurrence of distinct flutter bending modes, which 
differ from the ordinary vibration modes in that they 
tend to assume a form which approaches the next 
higher vibration mode and exhibit a correspondingly 
higher frequency. The flutter speed is consequently 



lower than that calculated on the basis of the lowest 
vibration frequency and the flutter frequency itself is 
higher. For ordinary damped structures, this effect 
lowers the flutter speed calculated on the basis of the 
lowest bending mode by only a few percent. 

3. A cantilever wing flutters at a speed calculated 
by using the constants for the most representative 
section, which is located at approximately three- 
quarters of the semispan. 

4. Aspect ratio and structural damping effects tend 
to increase the flutter speed by a few percent above 
that calculated for infinite aspect ratio and zero internal 
damping. 

5. Tin 4 effect of mass balancing to bring the center 
of gravity forward is essentially as predicted by theory. 
The effect of nacelles is of lesser importance, but large 
weights located at some distance away from the wing 
and attached to it show a very detrimental effect on 
the flutter speed. 

6. Wing-aileron experimental studies show that the 
characteristic flutter range predicted by the theory 
exists and is in substantial agreement with the pre- 
dicted values. A decrease in the unbalance and an 
increase in the frequency ratio are both beneficial. 
There exists, for each value of the unbalance, a certain 
critical frequency ratio and, inversely, for each fre- 
quency ratio, a certain critical value of the unbalance. 

7. The considerable difficulty involved in the deter- 
mination of the primary structural parameters includ- 
ing the damping is recognized and will constitute one 
of tin 4 chief problems of future flutter research. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., September 22, 1938. 
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LIST OF NOTATION 

a, angle of attack (fig. 40). 
0, aileron angle (fig. 40). 

h, vortical distance (fig. 40). 
da •• d 2 a 

a 0 , ft,, A 0 , <p\, (f2, amplitudes and phase angles. 

b, half chord, used as reference unit length. 

a, coordinate of axis of rotation (torsional axis) 
(fig. 40). Location of stiffness axis in percentage 
total chord measured from the leading edge is 



100^ or 



2 (stiffness axis) 



-1. 



2 w " 100 
coordinate of aileron hinge axis (fig. 40). Location 

total chord 
1+c 



•1. 



p, 
M, 



of aileron hinge axis in percentage 

measured from leading edge is 100— 

_ 2 (aileron hinge) 
C ~ 100 
mass of air per unit of volume, 
mass of wing per unit span length. 

7 2 

-r^p the ratio of the mass of a cylinder of air of a 

diameter equal to the chord of the wing to 
the mass of the wing, both taken for equal 
length along the span; this ratio may be 
expressed as k=0.24(6 2 /TT t ) (p/p 0 ) where Wis 
weight in pounds per foot span, b is in feet, 
and p/po is ratio of air density to standard 
air. [The quantity tr=TP/4& a (weight per 
square foot per chord in feet) has been used 
by British writers. Thus, k= (0.06/<t) (p/p 0 )J 

location of center of gravity of wing-aileron 
system measured from a (fig. 40); S a , static 
moment of wing-aileron per unit span length 
referred to a. Location of center of gravity 
in percentage total chord measured from the 
l+a+x a 



Mb' 



leading edge is 100- 

2 (center of gravity) 
100 



or a-\-x a - 



= Alb' rec ^ uce( ^ l° ca ^i° n °f center of gravity of aileron 
referred to c (fig. 40). S^, static moment of 
aileron per unit span length referred to c. 3 



"\Mb 2 



yi radius of gyration of whig aileron referred 
to a (fig. 40). I a , moment of inertia of 
wing aileron about the elastic axis per unit 
span length. 

^Jg, rrf « rfrail i,,sof gy r« tiM ,,,f,,il,™,re t .rre (1 
to c (fig. 40). 7/3, moment of inertia of 
aileron about c per unit span length. 

3 Note that M refers to the total wing mass and not to the mass of the aileron alone. 
40 



C a , torsional stiffness of wing around a per unit span 
length. 

Cp, torsional stiffness of aileron around c per unit 
span length. 

C h , stiffness of wing in deflection per unit span lengt h. 

Wa= vT' natural angular frequency of torsional vibra- 
tion around a in vacuum (u a =2?r/a, where 
f a is in cycles per second). 

V® 
-j) natural angular frequency of torsional vibra- 
0 tions of aileron about c. 



00 0 



Wft= -*/^™ natural angular frequency of wing in deflec- 
tion. 



Lead "ma Quart er Trailing 
edae chard M,d chord edg * 

ri 

/ /Axis of rota tiop 




e.g. of oi/cran - - - 



Figure 40.— Half chord b is used as the unit length. The positive directions of a, 
0, and h are indicated by arrows. Note that a is measured from raidchord and 
that Xa is measured from the elastic axis positive to the right. Also note that rp 
is a "reduced" parameter and not the actual distance from the hinge to the c. g. of 
the aileron. 

t 1 time. 

v, speed of forward motion. 
Vf, flutter or critical speed. 
o), circular frequency of wing vibrations. 

fc=— ) reduced frequency = number of waves in the 

wake in a distance equal to the semichordX27r. 

l/k, reduced wave length = length of one wave of the 
wake in terms of a distance equal to the semi- 
chord X2?r. 



F and G, functions of k in table 2. 



\2F 



R ah = —A h i+(^+ 



R»„ = -£«,- |^ 2 [(i - off- jfj 

n I 1 D T l2 /T n 2G T\,2F\ 





-B hX 




2G 
k 


Rc<x = 


— Gal 


-a 


-")t + w 


Rc&= 




2t 


2G,T 10 2F 
k tt k ' 2 


Rch= 


— Chi 


2G 
k 





I aa 



i{-(^»)f-a-«'>^ ■-■-] 
'•.-i[-G+«> / i 

'•.=r[i(v'f+f; 2f )+ B «] 
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v/bco h , flutter-speed coefficient (case 2). 

coa/W, frequency ratio (case 1). 

cop/wa, frequency ratio (case 2). 

cop/coa, frequency ratio (case 3). 

"7 (i/9)+a divergence velocity. 

gr a , #0, <7„, structural damping coefficients; wg corre- 
sponds approximately to the usual logarithmic 
decrement. 

M„, Mp, magnitude of sinusoidal impressed torques in 

the a and ft degrees of freedom. 
P 0 . magnitude of impressed force in the h degree of 

freedom. 
\p, phase angle. 

«b/««i, Po/Pst, ho/hst, peak response for the various 
degrees of freedom. 
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K \ CO / 



lb 0 = Ibd + fffi - h " * 
I ch — I eh + Ql$hX 

The quantities A al , Ap u etc. and Q a X, QpX, etc. are 
defined under the calculation scheme (pp. 5-7). The 
Ts are listed in table 2. The definitions of the T's 
are given in reference 1, page 5, should other values 
than those listed in the table be required. 
v/booa, flutter-speed coefficient (cases 1 and 3). 



TABLE 1 



L0 

6 
4 
3 
2 

1.5 

1.2 

1.0 
.80 
.66 
.60 
.56 
.60 
.44 
.40 
.34 
.30 
.24 
.20 
.16 
.12 
. 10 
.08 
.06 
.05 
.04 
.025 
.01 

0 



1/k 



0. (KK) 
. 100 

. Ui607 
.250 
. 33333 
.500 
. 66667 
. 83333 
1.000 
1.250 
1.51516 
L. fifi667 

1 . 78572 
2.000 

2. J7J7:i 
2. 500 

2. I'll I s 
3. 33333 
(. 16667 
5.000 
6. 250 
8. 33333 
10.000 

12. : 

16. 66667 
20. 000 
25.000 
40.000 
100.000 



I). .")( M III 
. -.006 

.5017 
.5037 
. 5063 
.5129 
.5210 
.5300 
.5394 
. 55 1 1 

. 5699 
.5788 
. 5S57 

. 5979 
. 6130 
. 6250 
. 6160 
.6650 
.60S! I 
. 7276 
.7628 
.8063 
.8320 

. S601 

. *020 
. 9090 
. 9267 
.9545 
.9S21 
1.00(1 



.0124 
. 0206 
.0305 
.0400 
.0577 
. 0736 
. OS77 
. 1003 
. 1165 

. ims 

.1378 
. I 12S 
.1507 
. 1592 
. 1650 
. 1 

. 1793 
. 1S62 
. 1 SSI] 
. 1876 
. LS01 
. 1723 
. 1604 
. I 126 
. 1305 
. 1160 
.0872 
.0482 



-2G/k 



0 

.00248 

.00686 

. (11525 

.02667 

.0577 

. 00 Is 

. 1 102 

. 2006 

.2912 

.3964 

.4593 

.5100 

.6028 

.7236 

. slV,0 

1.022 

L 195 

1.552 

1.886 

2. 345 

3. 002 
116 

4.010 

4. 753 

5. 220 
5. 800 
6. 976 
9. 640 



2F/& 



0 

.010012 
. 02787 
. 06296 
. 1125 
. 2565 
.4631 
.7361 
1.0788 
1.7316 
2.6166 
3. 2156 
3. 7353 
1. 7S32 

6. 3326 

7. 8125 
11. 192 
14. 778 
21.267 
36. 380 
50. 502 

111.99 
166. 4 
268.9 
105. 6 
727.2 
1158.3 
3051. 1 
19648 
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TABLE 2 
Values of T 



-1 


-0.5 


0 


0.1 


0.2 


0.3 


-1.00000 
-1.25000 
-1.00000 


-0. 54008 
-.38882 
-.80450 


-0. 21221 
-.08191 
-.50000 


-0. 16539 
-.05434 
-. 43644 


-0. 12490 
-.03429 
-.37353 


-0. 09064 
-. 02031 
-.31192 


-1.12500 
1.00000 
I. 50000 
.50000 

0 


-. 37922 
. 94233 
1.01125 
. 20075 

-.06891 


-.06250 
.81831 
.5os3i 

.00831 
-. 10610 


-.03540 
. 7Kis3 
.48813 
.05168 

-. 10452 


-.01672 
. 74779 
.41140 
. 03793 

-.09980 


-.00489 
. 70666 
.33870 
. 02679 

-.09105 


0 


. 13783 


.31831 


. 34789 


. 37426 


. 39474 


1.50000 


.87342 


.35611 


.27909 


.21187 


. 15555 


1.50000 


. 81355 


.28416 


.21304 


. 15369 


. 10565 


0 


. 05376 


. 05783 


. 05275 


. 04643 


. 03922 



T,/x 

TzI-k* 

T<M= -Cm).. 

Tih 

T,o/*- 

T„/2ir 

Tu/2k 

p/» — 

Aft- — - — 

Bf»- —. -. 





0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


Tilx _ 


-0. 06238 


-0. 04008 


-0. 02322 


-0.01145 


-0.00420 


-0.00076 


0 


T 3 M 


-.01107 


-.00539 


-.00223 


-.00071 


-.00014 


-.00001 


0 


Tilx( = -C fi 2) ---- 


-. 25231 


-. 19550 


-. 14238 


-.09400 


-.05204 


-.01808 


0 


Tifr 


-.00161 


. 00421 


.00428 


. 00302 


.00144 


.00032 


0 


TV*-.... 


. 06075 


. 00900 


. 549*2 


.48050 


. 39582 


.28231 


0 


Tn\2x 


.27029 


. 20075 


. 14874 


.09712 


. 05314 


.01889 


0 


T l2 /2x 


. 01798 


.01125 


. 00030 


. 00300 


.00110 


. 00020 


0 


Pi* 


-.08169 


-.06891 


-.05433 


-.03864 


-.02292 


-.00879 


0 


^ 


. 40844 


.41350 


. 40744 


.38644 


. 34378 


.26843 


0 




.10685 


.06892 


.04008 


. 01984 


.00830 


.00131 


0 




.06820 


. 04042 


.02118 


.00914 


.00277 


.00035 


0 


B„- T >-™ 


















. 03156 


.02384 


. 01651 


.00999 


.00476 


. 00127 


0 



The expressions for the T's arc listed in reference 1, page 5. 



TABLE 3 
Values of R aa " 



Raa — A a \ "4" Ran" 





-0.5 


-0.45 


-0.4 


-0.3 


-0.2 


-0.1 


0 


0 


0 


0 


0 


0 


0 


0 


0 


. 100 


0 


-. 00001840 


-.00122 


-. 00240 


-. 00352 


00400 


-. 00563 


. 16607 


0 


-. 0017193 


-. 0034 1 


-.00567 


-. 00980 


-. 01079 


-. 01566 


.250 


0 


-. 0038724 


-. 00707 


-. 01503 


-. 02209 


-. 02884 


-. 03539 


. 33333 


0 


-. 0068918 


-. 01305 


-. 02677 


-. 03935 


-.05140 


-. 06292 


.500 


0 


-. 015566 


-. 03084 


-. 06053 


-. 08907 


-.11645 


-. 14268 


. 66007 


0 


-. 027658 


-.05484 


-. 10779 


-. 15884 


-. 20799 


-. 25525 


. 83333 


0 


-. 043749 


-. 08677 


-. 17061 


-. 25153 


-. 32953 


-. 40460 


1.000 


0 


-. 063468 


12593 


-. 24786 


-. 36577 


-. 47966 


-. 58955 


I. 250 


0 


-. 10041 


-. 19937 


-. 39291 


-. 58003 


-. 76253 


-. 93860 


1.51516 


0 


-. 14966 


-. 29734 


-. 58674 


-. 86822 


-1. 14178 


-1.40740 


1.00607 


0 


-. 18260 


-. 36290 


-. 71661 


-1.00113 


-1.39647 


-1. 72263 


1.78572 


0 


-. 21098 


-. 41943 


-. 82866 


-1. 22769 


-1.61052 


-1.99515 


2. 000 


0 


-. 26779 


-. 53257 


-1.0.5309 


-1.56125 


-2.05795 


-2. 54230 


2. 27273 


0 


-. 35100 


-. 69838 


-1.38230 


-2. 05174 


-2. 70670 


-3. 34720 


2.500 


0 


-. 42981 


-. 85555 


-1.69450 


-2. 51700 


-3. 32300 


-4. 11250 


2.94118 


0 


-. 60814 


-1.21120 


-2. 40192 


-3. 57222 


-4. 72208 


-5. 85150 


3. 33333 


0 


-. 79560 


-1.58.535 


-3. 14680 


-4. 68435 


-6. 198(H) 


-7. 68775 


4. 16667 


0 


-1.2871 


-2. 56638 


-5. 10072 


-7. 00002 


-10.07728 


-12. 52150 


5.000 


0 


-1.9086 


-3. 80774 


-7. 57776 


-11.31006 


-15.00404 


-18. 66150 


(5. 250 


0 


-3. 0910 


-6. 17025 


-12.293(50 


-18. 37005 


-24. 39960 


-30. 38225 


8. 33333 


0 


-5. 7421 


-11.46918 


-22. 87832 


-34. 22742 


-45. 51648 


-56. 745.50 


10.000 


0 


-8. 4837 


-16. 95014 


-33. 83136 


-50. 64366 


-67. 38704 


-84. 06150 


12.500 


0 


-13. 635 


-27. 25090 


-54. 42160 


-81. 51210 


-108.52240 


-135.45250 


16. 66667 


0 


-25. 006 


-49. 98777 


-99. 88048 


-149. 67813 


-199. 38072 


-248. 98825 


20. 000 


0 


-36. 608 


-73. 18980 


-146. 27520 


-219. 25620 


-292. 13280 


-364. 90500 
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TABLE 3— Continued 

ValUCS Of laa" 



I aa 1 1 h I <ui 





0 5 


—0 45 


O 4 


Q y 


— 0.2 


q j 


0 


0 


1. 00000 


0. 90250 


0.81000 


0. 04000 


0. 49000 


0. 36000 


0. 25000 


. 100 


i. ooooo 


. 90256 


. 81014 


. 04031 


. 49044J 


. 36070 


. 25094 


. 16067 


1. 00000 


. 90208 


. 81038 


. 04083 


.49135 


. 30192 


. 25258 


. 250 


1. ooooo 


. 90291 


. 81086 


.04187 


. 49303 


. 36432 


. 25578 


. 33333 


L. 000*00 


. 90323 


. 81164 


. 64331 


. 49535 


. 36764 


. 26019 


.500 


1. ooooo 


. 90415 


. 81345 


.04741 


. 50189 


. 37689 


. 27240 


. 00007 


1.00000 


. 90524 


.81570 


. 05224 


. 50962 


. 38784 


. 28690 


. 8.3333 


l.(XXXK) 


. 90690 


.81922 


. 66002 


.52120 


. 40404 


.30810 


1 000 


1 1 II II II II 1 


in )s -'( 


82297 


66751 


. 53363 


. 42133 


. 33000 


f. 250 


L ooooo 


! 91192 


.' 82938 


! 68003 


. 55464 


.45051 


. 30855 


1.51510 


1.00000 


.91568 


. 83706 


.09091 


. 57956 


. 48501 


. 40325 


1.00007 


I. 001)00 


. 91797 


.84175 


. 70664 


. 59469 


. 50590 


. 44025 


I. 78572 


1. ooooo 


.91986 


. 84557 


.71457 


.60701 


. 52280 


.46215 


2. 000 


1.0000(1 


. 9233 i 


. 85266. 


. 72923 


. 62972 


.55413 


. 50245 


2. 27273 


1.00000 


. 92794 


. 80202 


. 74856 


. 65062 


. 59520 


. 55530 


2. 500 


1.00000 


. 93187 


. 87000 


. 76500 


. 68500 


. 63000 


. 60000 


2. 94118 


1. ooooo 


. 93964 


. 88570 


. 79739 


.73190 


. 09829 


. 68755 


3. 33333 


1. 00(X/0 


. 94657 


. 89980 


. 82620 


. 77920 


. 75960 


.76,500 


4. 10007 


1.00000 


.96120 


. 92949 


. 88675 


. 87200 


. 88533 


. 92655 


5. 000 


1. ooooo 


. 97518 


. 95763 


. 94437 


.90021 


1.00515 


1.07D20 


6. 250 


1.00000 


. 99478 


. 99720 


1.02490 


1.08312 


1. 17180 


L 29110 


8. 33333 


1.00000 


1.0235 


1. 05507 


1.24238 


1.26195 


1.41378 


L 59785 


10. 000 


1.00000 


1. 0433 


1.09484 


1. 32296 


1. 38436 


L 57904 


1. 80700 


12.500 


1.00000 


1.0688 


1. 14613 


1.42667 


1 54163 


1.79101 


2. 07480 


10. 60007 


1.00000 


I 1029 


1.21474 


1. 56516 


1.75126 


2. 07304 


2. 43050 


20. U00 


L 00000 


L 1240 


1. 25838 


1.65312 


1. 8S422 


2. 25168 


2. 05550 



TABLE 4 
Values of Rbp" 
R^-Bh + Rv" 





0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.0 


0.7 


0.8 


0.9 


1.0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


.100 


.00123 


.00100 


.00078 


. 00060 


. 00045 


. 00032 


. 00021 


.00011 


. 00005 


.00001 


0 


.16667 


. 00342 


.00207 


. 00219 


.00168 


. 00124 


.00087 


. 00057 


. 00032 


.00014 


. 00004 


0 


.250 


. 00768 


. 00623 


. 00493 


. 00378 


.00279 


. 00190 


.00126 


.00071 


. 00032 


. 00008 


0 


. 33333 


.01367 


.01 LOS 


.OOS77 


. 00673 


. 00498 


.003 lit 


.00225 


.00119 


. 00056 


.00015 


0 


.500 


. 03085 


. 02505 


. 01978 


.01519 


.01122 


.007S7 


. 00507 


. 00290 


.00126 


. 00033 


0 


. 66667 


. 05494 


.04401 


. 03525 


. 02700 


.01999 


. 01402 


.00005 


.00515 


. 00225 


. 00059 


0 


. 83333 


.08647 


.07018 


. 05540 


.04250 


. 03137 


.02230 


.01419 


.OOM 10 


. 00352 


.00002 


0 


1.000 


. 12517 


. 10158 


. 08015 


. 06146 


.04535 


.03177 


.ii2i Ms 


.01104 


. 00510 


.00133 


0 


1.250 


. 19729 


. 16000 


. 12021 


. 09670 


. 07130 


.04990 


. 03214 


.01S25 


. 00798 


.0020S 


0 


1.51510 


.29270 


. 23723 


. 18099 


. 14317 


. 10547 


. 07373 


. 04744 


. 02690 


.01175 


. 00307 


0 


1. 66667 


. 35014 


. 28854 


. 22733 


. 17398 


.12810 


.08952 


. 05755 


. 03262 


. 01425 


.00371 


0 


1. 78572 


. 410G2 


. 3325S 


. 26195 


. 21 id 111 


. 14750 


. 10302 


.00020 


.03750 


.01637 


. 00426 


0 


2. 000 


.51910 


. 12H22 


. 33078 


.25290 


. 18595 


. 12980 


.08335 


.04717 


.0205s 


. 00535 


0 


2. 27273 


. 67692 


. 54753 


.43071 


. 32902 


.24177 


. 10S57 


. 10813 


.06113 


.02001 


.00601 


0 


2. 500 


. 82550 


. 00738 


. 52403 


. 10050 


. 2U407 


. 20488 


. 13131 


.07417 


. 03230 


. 00838 


0 


2.94118 


1. 151)07 


. 93605 


. 73500 


. 56040 


.41094 


.28588 


. 18290 


. 10318 


. 04486 


.01162 


0 


3. 33333 


1.50670 


1.21566 


. 95304 


. 72637 


.53204 


. 36968 


. 23630 


. 13308 


. 05779 


.01495 


0 


4. 10667 


2. 40775 


1.93923 


1. 51852 


L 15435 


. 84375 


.58495 


. 37303 


.2095S 


.09080 


. 02342 


0 


5.000 


3. 53385 


2. 84190 


2. 22193 


I.ON02S 


1. 23036 


.85134 


.54187 


. 30380 


. 13137 


.03380 


0 


6.250 


5. 65191 


4. 00926 


3. 54031 


2. 08138 


1.95216 


1. 34764 


. 85572 


.47853 


. 20644 


.05296 


0 


8. 33333 


10. 34067 


8.28124 


6. 44720 


4. 87081 


3. 53070 


2. 43454 


1.54135 


. 85928 


. 36960 


. 09447 


0 


10. 000 


15. 14353 


12. 11112 


9. 41589 


7. 10321 


5. 14958 


3. 53887 


2. 23667 


L 24464 


. 53443 


. 13633 


0 



TABLE 4— Continued 
Values of I b &" 



l/*\^ 


0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


0 


0. 32298 


0. 23827 


0. 16930 


0.11472 


0. 0730G 


0. 04275 


0. 02213 


0. 00944 


0.002S3 


0. 00035 


0 


.100 


. 322S9 


.23820 


. 16925 


.11408 


.07304 


. 04273 


.02212 


.00044 


.00283 


. 00035 


0 


. 10(567 


. 32273 


. 23808 


. 16916 


. 1 1402 


. 07300 


. 04269 


.02211 


. 00943 


. 00283 


. 00035 


0 


.250 


. 32212 


. 23784 


. 10S09 


.11450 


.07292 


.04200 


. 02209 


. 00942 


.00282 


. 00035 


o 


. 33333 


. 32198 


. 23751 


. 16874 


. 11433 


.07280 


. 04243 


. 02205 


.00940 


.00282 


.00035 


o 


.500 


.32075 


.23658 


. 10SOO 


.11386 


. 07250 


.04241 


.02195 


. 00937 


. 002SO 


. 00035 


0 


. 66667 


.31931 


. 23548 


. 16727 


.11331 


. 07213 


. 04220 


.02184 


.00931 


. 00279 


. 00034 


0 


. 83333 


.31714 


. 233.S5 


. 16609 


.11249 


.07161 


. 04183 


.02108 


. 00925 


. 00277 


. 00034 


0 


1.000 


.31483 


.23212 


. 16484 


. 11163 


. 07106 


.04156 


.02150 


. 00917 


.00274 


. 00034 


0 


1.250 


. 31090 


.22919 


. 16273 


.11019 


. 07013 


.04101 


.02121 


. 00904 


.00271 


. 00034 


0 


1.51516 


. 30625 


. 22572 


. 16024 


. 10849 


. 06903 


. 04036 


.02087 


. 00890 


. 00267 


. 00033 


0 


1. 66667 


. 30343 


. 22362 


. 15S73 


. 10710 


.06837 


. 03997 


. 02067 


.00881 


. 00264 


. 00033 


0 


1.78572 


.30112 


.22190 


. 15752 


. 10662 


.00783 


. 03966 


.(12(150 


. 00874 


. 00262 


. 00033 


0 


2. 000 


. 290ss 


.21870 


. 15526 


.10509 


. 06685 


.03908 


.02021 


. 00861 


. 00257 


. 00033 


0 


2. 27273 


.29130 


. 21462 


. 15231 


. 10307 


.06556 


. 03832 


.01981 


. 00845 


. 00252 


. 00032 


0 


2.500 


.28657 


.21112 


. 1 10SO 


.10137 


.00148 


.03708 


. 01948 


. 00831 


. 00249 


.00031 


0 


2. 94118 


. 27726 


. 20423 


. 14491 


.OOM 13 


. 06235 


. 03043 


.01883 


.00803 


. 00240 


. 00030 


0 


3. 33333 


. 26899 


. 19813 


. 14056 


.0950!) 


. 06046 


. 03533 


. 01826 


. 00778 


. 00232 


. 00029 


0 


4. 16667 


. 25166 


. 18536 


.13152 


.08891 


. 05655 


. 03305 


. 01708 


.00728 


.00217 


. 00027 


0 


5.000 


. 23522 


. 17325 


. 12292 


.08315 


.05286 


.03089 


. 01596 


.00681 


.00203 


. 00025 


0 


6. 250 


. 21229 


.15642 


.11100 


.07510 


. 04775 


.02791 


.01442 


.00615 


.00183 


.00023 


0 


8. 33333 


. 17895 


.13197 


. 09372 


.06345 


.04038 


. 02362 


. 01220 


. 00521 


.00155 


.00019 


0 


10. 000 


.15613 


. 1 1525 


.08194 


.05551 


. 03535 


. 02069 


. 01070 


.00457 


.00135 


.00017 


0 
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2 



1.09 



6.8 



19.7 



15.0 



1.29 



21.0 



16.7 
17.0 

14.8 



2B 



1.31 



17.8 



c. g. location 

(percent 
chord from 
L. E.) 



12. .-> 



11.2 



42.3 



9 2 



2 c 



0.3125 



0. 3125 



. 32 



30. 50 



35. 95 



Test conditions and remarks 



Tests on wing 1 were preliminary. 

A stop w as used. 
Later tests employed a screen instead 

Damping test.. . * 

Run 20-3 repeated... 

Weight, 0.028 lb. on L. E. 18 in. from 

tip. 



Stop used; all other tests screen used 

(tip nutter). 

Model support fractured 

Model repaired 

Restraining wires 17 in. from tip, 

3.6 in. from L. E. 
Nacelle 6.10 lb. placed 2 ft. from top, 

2.70 in. ahead of c. g. 
Nacelle placed 2 ft. from top, 1.34 in. 

back of c. g. 

Run 20-16 repeated 

Counterweight 3.07 lb. at L. E. 13 in. 

from tip. 

Run 20-20 repeated 



Plain; no weight attached 

Disk 12.77 lb. 12-in. diam. on bottom; 
c. g. of disk on c. g. of airfoil. 

Disk 12.45 lb. 17-in. diam. in same po- 
sition as 20-24. 

Disk as in 20-25 placed 1 in. ahead of 
c. g. 

Disk as in 20-25 placed 1 in. back of 
c. g. 

Run 20-23 repeated. 

Airfoil tip restrained by wires 

Run 20-25 repeated _ 

Run 20-23 repeated 

Same as 20-31, movies taken 

Plain; restraining wires 8 in. from tip. 
Plain; restraining wires 60 in. from tip. 
Restraining wires 40 in. from tip; 
model broke near tip. 



20-2 

20-3 
20-4 
20-5 
20-8 



20-11 
20 15 
20-16 

20-17 

20-18 

20-19 
20-20 



20-23 

20- 21 

20-25 

20-26 

20-27 

20-28 
20-29 
20-30 
20-31 
20-32 
20-33 
20-34 
20-35 



2.283X10-3 
. 278 



2.274 
2. 265 



2. 154 
2.216 



2.116 
2.170 



2. 225 
2. 173 



2.174 



2. 159 
2.204 

2. 228 

2.232 

2. 214 

2. 149 
2. 154 
2. 232 
2.146 
2.144 
2. 122 
2.171 
2. 053 



15.0 



10.2 
9. 45 



10.2 
8. 35 

7. 35 

7. 48 

7. 25 



10.8 
13.1 
8. 24 



MECHANISM OF FLUTTER 45 



TABLE I.— EXPERIMENTAL RESULTS OF FLUTTER INVESTIGATION — Continued 



Figure number 


Wing model 


Wing dimensions 


First bending (cycles/sec.) 


Second bending (cycles/sec.) 


Third bending (cycles/sec.) 


Fourth bending (cycles/sec.) 


First torsion (cycles/sec.) 


Second toision (cycles/sec.) 


c. g. location 

(percent 
chord from 
L. E.) 


Stiffness axis (peicent chord 
from L. E.) 


& 
.2 

29 

j~ o 

O o 

w a 
9 

1 

s 


Wing weight (lb.) 


Test conditions and remarks 


d 
fA 


Air density, p 


Flutter frequency (cycles/ 
sec.) 


Flutter velocity (m. p. h.) 


Theoretical results 


O 
09 

1 

o 


to 

> 

i 

8 

o 


11 


2C 


Rectangular: 12-in. 
chord; 6-ft. 4-in. 
span 


1.40 
1.22 

1.29 
1.34 
1.38 
1.39 
1.39 
1.21 


8.2 
7.2 

8.2 
8.0 
7.6 
7.8 
8.0 
7. 7 


22.5 
20.8 

22.4 
22.2 
22.5 
22.0 
21.2 
19.8 




18.8 
16.2 

15.8 
16.0 
16.8 
17.4 
18.2 
13. 8 




42.5 




30.0 







Model 2B repaired 


20-73 
20-74 

20-75 
20 76 
2(1 77 
20-78 
20 7'.i 
20-81 


2.133XHH 
2. 169 

2.152 
2. 101 
2. 119 
2. 122 
2. 136 
2. 056 




220.4 

MM'., s 

210. 5 
244. 6 
235. 7 
232.2 
2 Ml. s 
294.3 


CO 












With 3.u5-lb. counterweight 1 in. from 

tip at L. E. 
Counterweight 13 in. from tip at L. E . 
( 'oiinteru eight 25 in. from t ip at t>. K 
Counterweight 37 in. from tip at L. E. 
( 'oimterweighi 1!) in. from tip at L. V. 
Counterweight 61 in. from tip at L. E. 
With weights simultaneously at 13 in.. 

37 in., and 61 in. from tip. 




















































































































12 


3 
4 


CN 
O 

oa 
>-i 
u 

& 

03 


1. 65 
1.64 


6. 8 
5.2 


16. 5 
13.1 


30. 7 
28.4 


29. 1 
8.9 

O 1 

y. i 
8.8 





42.5 


42.0 


30.0 




20.48 


Plain tapered wing 


20-22 
20 68 
20-69 

20-70 

20-71 
20-72 


2.160 
2. 150 
2. 232 

2. 245 

2.240 
2. 149 


16.4 


216.0 
214.9 
109.2 

91. 0 

90.4 
213.3 




37.8 




Run 2U-22 repeated . 


Nacelle 6.10 lb.; c. g. 3.10 in. back of 

c. g. of airfoil. 
Nacelle placed 1.10 in. back of c. g. 

Fig. 25. 

Nacelle placed 2.80 in. ahead of c. g-.. 
Ruo 20-68 repeated.- 




























8.50 
14.0 














13 


2. 07 


6.0 


13.3 


24.0 


36.2 




41.2 


41.5 






16.88 


Violent flutter suddenly. Wing bent 
out of shape. Fig. 20. 


20-67 


2.174 




203.9 












12 


3 


g 

u 

& 

a 
H 


1.65 


6.8 


16.5 


30.7 


29.1 




42.5 


42.0 






20. 48 


Restraining wires 59 in. from tip at 
stiffness axis. 

Same except wires 2 in. ahead of stiff- 
ness axis. 

Same except w ires 1 in. back of si ill- 
ness axis. 

Restraining wires 40 in. from tip at 

stiffness axis 
Restraining wires 17 in. from tip at 

stiffness axis. 
Restraining wires 1 in. from tip at 

stiffness axis. 
Restraining wires 27 in. from tip at 

stiffness axis. 


20 Ml 
20-90 
20-91 
20-93 
20-94 
20-95 
20-96 


2. 187 
2.200 




211.1 
204.4 


CS 

1 ti 

£ 










2. 234 
2.229 
2.182 
2. 202 




155.8 
152.9 
211.4 
200.8 


16 


6 




2.30 


12.8 


31.2 




19.2 




42.5 


45.2 






3. 45 


Rectangular wooden model; 81-in. 
span. 


20-82 


2. 268 




73.8 












16 


7 




4.92 


25.0 


67.2 




35.8 




42.5 


43.9 






6. 01 


Tapered wooden model; completely 
destroyed by violent flutter. 


20-92 


2.199 




205.6 
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Wing 
model 



Aileron (see fig. 14) 



Wing 
weight 
(lb.) 



Aileron 
weight 
(lb.) 



First 
bend- 
ing 

(cy- 
cles/ 
sec.) 



Aileron 

fre- 
quence 
(cy- 
cles/ 
sec.) 



r fi * (cf. 

nuta- 
tion) 



Test condit ions and remarks 



Run 



Air dens- 
ity, p 



Lower 
critical 
velocity 
(m. p. 
h.) 



Air dens- 
ity, p 



AIL. 
AV. 
AIL 
AV 

Ml 
All 
AIL. 

Ml 

Ml 

Ml 
All 

Ml 

All 

All 

All 

AIL 

All 



All 
All 

AIL— 
MI 
All 

AI 

A V 1 1 
.Mil 
All J 

Am... 

AIIL.. 
\ 1 1 1 
All! 
AVI I I 
AVIIL 
Will 
AVI I I 
.Will 



44.84 



n. '.mi 
.618 

.05(1 

.nis 
.950 



10. 68 



.567 
.300 
1.43 



2.01 



0 
0 

5.75 
7.91 
12.5 
s. :»() 
11.0 
11.0 
12.5 
13.1 
9. 67 
9.17 
5.75 
0 

13.3 
10.8 
0 

0 
0 
0 

3.00 

6. 92 

0 

0 

0 

0 

0 

0 

4.67 

6.00 

0 

0 

0 

0 

7.00 



(I. nun. 



Aileron with 3 free hinges. 



Aileron with springs - 

Theoretical results for aileron All shown in fig. 26 



Hinges free.. 



3 free hinges with aileron counterweight, 0.557 
lb., Fig. 15. 

Counterweight, 0.452 11).... -- - . 

Counterweight, 0.346 lb... 

Counterweight, 0.398 lb 

.....do. 



do 

2 free hinges 

3 free hinges. ._ 

4 free hinges 

4 free hinges with lead for damping . 

do. 

do 

do 



.do. 



Hinge-pin i>osition varied as shown in fig. 27. 

do _ - — 

do - 

.—do 

do 



20-36 
20-37 
20-38 
20-39 
20-40 
20-41 
20 12 
20-43 
20 11 
20-45 
20-46 
20-47 
20-48 
20-49 
20-50 
20-51 
20-52 

20-63 
20-54 
20-55 
20-56 
20-57 
20-58 
20-59 
20-60 
20-61 
20-63 
20-64 
20-65 
20-66 
20 s:\ 
20-84 
20 SO 
20-87 

20 hS 



2. 273X10-3 
2. 278 
2. 273 
2.271 
2.269 
2. 306 
2.301 



2. 259 
2. 2»3 

2. 245 
2. 242 
2. 242 
2. 240 



2. 249 
2.243 
2. 243 
2. 262 
2.262 
2. 250 
2. 250 
2. 251 
2. 261 
2. 255 
2. 250 
2. 260 
2.285 
2.281 

2. 250 
2.259 



53.9 
50.7 
50.4 
42. ■} 
76.8 
44.8 
55.4 



ov t 
57.4 
50.0 
52.6 
58.3 
60.3 
65.0 



2.240X10-* 
2. 229 



2. 236 
2. 162 
2. 165 
2. 138 
2. 200 
2.199 



No flutter 
No flutter 



75.2 
99.5 
99.9 
62.8 
65.3 
57.8 
60.2 
52.2 
78.0 
97.8 
67.3 
63.8 
64.4 
93.6 



2. 243 
2. 173 
2. 193 
2.218 
2.163 
2. 221 
2. 161 
2. 199 
2. 171 
2. 150 
2.178 
2.210 
2.178 
2.143 
No flutter 
85.5 I 2.098 
41.1 2.130 



290.8 
264.1 



C 




z 

Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 




Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Force 
(parallel 
to axis) 
symbol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal. 


X 


X 
Y 


Rolling 

Pitching 

Yawing 


L 


Y >Z 


Roll 


e 


ti 


V 


Lateral. _ 


Y 


M 


z — >x 


Pitch 


V 


Q 


Normal 


Z 


Z 


N 


X >Y 


Yaw 




w 


r 





















Absolute coefficients of moment 

V '~qbS Vm ~qcS 
(rolling) (pitching) 



(yawing) 



Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 



D, Diameter 

Geometric pitch 
pfD, Pitch ratio 
V, Inflow velocity 
V s , Slipstream velocity 

T 

T, Thrust, absolute coefficient C T = — ^th 
' pw u 

Q 



4. PROPELLER SYMBOLS 

P, 

c 9 , 



n, 



Power, absolute coefficient pn *jy> 

V~pV~ b 

Efficiency 

Revolutions per second, r.p.s. 
Effective helix an gl e=tan ~ 1 ^2^rrn) 



Q 9 Torque, absolute coefficient C Q = 



1 hp. = 76.04 kg-m/s=550 ft-lb./sec. 

1 metric horsepower— 1.0132 hp. 

1 m.p.h. = 0.4470 m.p.s. 

1 m.p.s. = 2.2369 m.p.h. 



'pn'D 5 

5. NUMERICAL RELATIONS 



1 lb. -0.4536 kg. 

1 kg=2.2046 lb. 

1 mi. = 1,609.35 m=5,280 ft. 

1 m=3.2808 ft. 



i 



HI 



o 



a? 

CO 

CO 



ZD 



j,4 V— (XL 

S ^ 



